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accelerometer 

Natural  frequency  of  sensor  ra 

Duaay  variable  used  in  linear 
analysis  _ 

Angular  rates  about  x and  y axes, 
respectively 


Units 


xv  i 


Slop*  of  Corqu*-sp**d  curve 


Slop*  of  torque-error  signal  curve  kg-a/rad 
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SECTION  1 

INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 

The  "Closed-Loop  Optimization  Program  for  the  M60A1  Tank  Gun 
Stabilization  System"  study  project  was  undertaken  by  the  Bendlx  Research 
Laboratories  under  Contract  DAAA09-74-C-2068  for  the  Rock  Island  Arsenal. 
The  objective  was  to  develop  a gun  stabilization  system  in  which  per- 
formance and  reliability  are  optimized.  It  had  been  concluded  prior  to 
the  initiation  of  this  contract  that  fluidic  approaches  appeared  attrac- 
tive from  the  reliability  and  cost  standpoint.  The  purpose  of  this 
contract  is  to  provide  the  Contracting  Officers  with  conclusive  evidence 
and  recommendations  regarding  the  applicability  of  the  fluidic  approach 
to  tank  weapons  systems,  while  supplying  a mathematical  model  of  the 
entire  vehide/control  system.  In  order  to  provide  this  information, 
fluidic  devices  were  to  be  considered  without  forcing  incorporation  of 
a particular  existing  fluidic  device. 

The  results  of  Phase  I,  discussed  in  this  report,  consist  of  a 
requirement  analysis,  a model  definition,  and  a system  stability  and 
performance  analysis.  A detailed  performance  study  of  the  different 
applicable  sensors  consisting  of  electronic  rate  gyros,  a hydraulic 
rate  sensor,  an  Integrating  accelerometer,  a laminar  vortex  sensor, 
and  a pneumatic  accelerometer  were  also  included  for  future  trade 
studies.  A complete  error  analysis  for  two  of  the  sensors  is  recom- 
mended in  order  to  augment  a decision  on  the  applicability  of  fluidic 
approaches  in  this  area. 

Control  system  design  was  performed  in  accordance  with  specifica- 
tions which  call  for  a pointing  accuracy  of  1/2  mil  rms  diameter  circle 
without  a gunner  operating  the  controls.  A 3 dB  bandwidth  of  at  least 
15  Hz  is  also  desired.  It  is  required  that  the  system  hardware  be  out 
of  saturation  98  to  99  percent  of  the  time.  The  duration  of  an  engage- 
ment was  specified  to  be  15  to  20  s. 

1.2  SCOPE 

for  the  purpose  of  deriving  a realistic  mathematical  model  of  the 
stabilization  system,  a study  was  made  of  existing  literature  on  the 
M60A1  tank,  and  of  previous  system  studies  on  the  vehicle.  Upon  simpli- 
fying models  of  the  vehicle  in  the  literature,  where  appropriate,  and 
independently  devloping  realistic  models  of  the  servovalve,  actuators 
and  sensors,  a complete  control  system  model  was  derived.  Selection 
of  the  appropriate  actuator  system  hardware  to  be  used  was  carried  out 
in  close  cooperation  with  the  Contracting  Officer's  Representatives. 
Models  were  defined  separately  for  the  azimuth  (turret)  and  elevation 


(gun)  axes.  In  addition  to  the  vehicle  and  gun  aodels,  the  Mathematical 
modal  included  the  controller  and  compensation,  the  effects  of  nonlinear 
valve  flow,  and  static,  coulomb,  and  viscous  friction  effects.  The 
affect  of  hull  dynamics  on  the  system  was  also  evaluated.  A mathematical 
modal  of  coupling  between  axes  due  to  roll  motion  was  investigated  analy- 
tically with  respect  to  tne  approximate  effects  on  system  performance 
and  pointing  accuracy.  Control  laws  were  derived  analytically  In  accord- 
ance with  the  performance  specifications  set  forth  prior  tu  the  start 
of  the  contract.  A control  law  was  specified  for  both  rate  and  position 
command  lnj  ut  concepts. 

from  t'ie  model  definition,  an  analog  computer  program  was  defined 
and  programmed  for  each  control  concept.  Diagrams  and  descriptions  of 
all  programs  used  In  this  studvare  contained  In  Appendices  A and  B.  The 
lock  Island  Arsenal  digital  program  "HITPRO"  was  adapted  to  the  Bendlx 
Research  Laboratories  computer,  and  tna  analog  programs  were, at  that 
point,  modified  to  hybrid  programs  in  order  to  accept  Inputs  from  HITPRO. 
These  Inputs  consisted  of  vehicle  rates  corresponding  to  those  experi- 
enced by  an  M60  tank  while  traversing  the  Aberdeen  Proving  Ground  terrain. 
Additional  hybrid  functions  of  thase  programs  consisted  of  a time  delay 
property  of  the  laminar  vortex  sensor. 

Prior  to  a complete  hybrid  computer  analysis  of  the  system,  an 
analytical  systems  study  was  conducted  In  order  to  develop  a fuller  com- 
prehension of  the  overall  system  behavior.  The  study  also  aided  In  the 
development  of  digital  check  solutions  for  all  of  the  hybrid  computer 
analyses.  In  addition,  the  analytical  study  Includes  a sensor  error 
analysis  and  an  evaluation  of  coupling  of  axes  on  pointing  errors. 

The  hybrid  computer  simulation  results  consist  of  am  evaluation  of 
the  effects  of  nonlinear  valve  flow  and  hull  dynamics,  a comparison  of 
the  rate  and  position  control  concepts,  and  a detailed  evaluation  of 
each  of  the  five  sensors  studied.  The  electric  rate  gyros,  as  they  are 
presently  used  In  the  field,  were  evaluated  In  both  tbs  elevation  and 
azimuth  axes.  The  hydraulic  rate  sensor,  the  Integrating  accelerometer, 
and  the  laminar  vortex  sensor  were  evaluated  In  tha  elevation  axis,  and 
the  pneumatic  accelerometer  was  evaluated  In  the  azimuth  axis. 

The  conclusions  reached  on  the  basis  of  the  computer  analysis  are 
contained  In  Section  7.  Recommendations  as  to  future  studies  and  stabi- 
lization approaches  are  contained  In  Section  8. 

1.3  SUMMARY 

A mathematical  model  of  a suitable  stabilization  system  for  the 
M60A1  tank  main  gun  was  formulated  and  programmed.  The  modal  was  formu- 
lated so  as  to  Include  most  of  the  significant  nonlinearities  such  as 
nonlinear  valve  flow  and  hull  dynamics  due  to  gun  motion.  A hybrid 
computer  analysis  was  performed  to  determine  the  operating  characteris- 
tics of  the  stabilization  system,  to  evaluate  prospective  sensors  for 
sensing  gun  and/or  hull  rate,  and  to  determine  whether  a rate  or  a 
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position  command  control  concept  is  preferable  with  respect  to  speci- 
fied performance  criteria. 

The  analytical  study  revealed  that  both  the  rate  and  the  position 
control  concepts  required  a proportional  plus  integral  control  lav  in 
order  to  minimize  the  gun  tracking  error.  It  was  also  shown  that  the 
rate  and  the  position  concepts  are  equivalent  in  terms  of  nulling  out 
the  effects  of  hull  motions,  and  thus  in  terms  of  stabilising  the  gun 
after  the  target  is  in  the  sight.  A computer  analysis  which  followed 
verified  this  equivalence.  In  addition,  it  was  possible  to  show  that 
the  effect  of  hull  motions  on  the  system  can  be  minimized  by  either  a 
high  control  loop  gain  along  with  a lead-lag  compensation  network  or 
by  using  a hull  sensor  signal  in  the  control  law. 

The  extensive  computer  simulation  analysis  revealed  several  signi- 
ficant conclusions  in  the  areas  of  stabilization  control  philosophy  and 
sensor  applicability.  In  the  process  of  arriving  at  a full  computer 
model  of  the  system  for  sensor  evaluation,  it  was  found  that  the  effect 
of  hull  dynamics  on  tha  gun  was  negligible.  Nonlinear  valve  flow,  how- 
ever, was  found  to  have  a significant  Influence  on  system  performance. 

A linearized  flow  model  was  not  sufficiently  accurate  for  use  in  this 
study. 

It  was  found  that  all  five  of  the  sensors  s’-udiad  meet  the  per- 
formance criteria  set  forth  by  the  Contracting  Officer's  Uepresentatlve. 
In  addition,  this  study  indicates  that  these  criteria  can  be  met  by  using 
only  a gun  sensor.  If  verified  by  further  studies,  the  need  for  a cor- 
responding hull  sensor  may  be  eliminated. 

The  detailed  sensor  study  revealed  that  automatic  offset  and.  inte- 
grator drift  nulling  circuits  are  required  when  using  an  acceleration 
sensor.  A method  which  can  be  used  for  this  purpose  is  described  in 
Section  5.5.  In  addition,  it  was  determined  that  increasing  the  gain 
of  the  acceleration  sensor  will  decrease  the  sensor  offset  effects  and, 
hence,  the  drift  rate.  More  generally,  in  the  sensor  study  it  was  found 
that  sensor  gain  errors  have  a null  effect  on  the  tracking  error.  Also, 
a combination  of  sensor  deadband  and  gun  or  turret  friction  will  cause 
the  system  to  limit  cycle. 

In  order  to  compensate  for  sensor  phase  lag,  feedback  compensation 
was  found  to  be  required.  Forward  path  compensation  is  desirable  for 
obtaining  stability  with  higher  loop  gains  for  this  system. 
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SECTION  2 

MATHEMATICAL  MODEL  DESCRIPTION 


2.1  GENERAL 

This  section  contains  t<ia  mathematical  model  of  the  gun  stabiliza- 
tion system  ss  programmed  on  the  hybrid  computer.  This  model  consists 
of  the  vehicle  sudel,  servovalve  and  motor,  sensor  models,  and  the  con- 
troller. The  azimuth  and  elevation  systems  are  presented  separately. 

The  effects  of  coupling  are  not  included.  The  individual  system  com- 
ponents sre  described  in  Sections  2.2  through  2.4.  Also  presented  are 
the  controller  andels  for  the  following  configurations: 

(1)  Rate  control  with  two  rate  sensors 

(2)  Rate  control  with  one  acceleration  sensor  in  gun  axis 

(3)  Position  control  with  two  rate  sensors 

(4)  Position  control  with  one  acceleration  sensor  in  gun  axis 

A simplified  block  diagram  of  a rate  control  system  with  two  rate  sensors 
Is  shown  in  Figure  2-1. 

2.2  VEHICLE  MODELS 

The  vehicle  andels  consist  of  an  elevation  (gun)  axis  model  and 
an  ezisRith  (turret)  axis  model.  They  Include,  respectively,  the  gun 
and  turret  inzrtias,  effects  of  friction,  and  hull  suspension  effects. 

2.2.1  Elevation  Axis  Vehicle  Model 

The  elevation  axis  vehicle  model  shown  in  Figure  2-2  con- 
sists of  the  combined  gun  and  hull  dynamics.  For  the  gun  dynamics,  the 
torque  applied  to  the  gun  equals  the  actuator  torque  T^,  plus  the  dis- 
turbance torque  Tj  acting  on  the  gun,  minus  the  trunnion  friction  torque 
Tfg;  It  follows  that  the  gun  angular  acceleration  with  respect  to  iner- 
tial space  is  given  by: 


Vr 


(tA  + Td  - Tf,> 


The  friction  model  is  discussed  in  Section  2.2.3.  The  disturbance  torque 
T4  is  generated,  for  example,  by_a  gun  firing.  The  hull  dynamics  gener- 
ate a hull  angular  acceleration  6^  with  respect  to  the  tracks.  This 
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Figure  2-2 . Elevation  Axis  Vehicle  Model 
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acceleration  results  from  the  reaction  of  the  actuator  torque  and  the 
tun  firing  torque  Tgf.  The  acceleration  la  given  by: 


+ T 


D 

aa 


+ V 


where  the  last  two  terms  represent  the  hull  suspension. 

The  terrain  pitch  rate, 8^  added  to  the  relative  hull  rats 
8t,y  gives  the  total  hull  pitch  rate  6^  as  shown  below. 


eh“ 


6.  + 
hv 


ht 


The  angular  rate  of  the  gun  with  respect  to  the  hull  (relative  gun  rate) 
is  given  by: 


A list  of  the  vehicle  elevation  parameters  is  contained 

in  Table  2-1. 

Table  2-1.  Elevation  Vehicle  Parameters 


Parameter 

Value 

Unite 

D.e 

7.8  x 104 

kg-m-s 

Dgh 

166 

kg-a-s 

J* 

527 

. 2 
kg-m-s 

Jhe 

1.73  x 104 

. 2 
kg-m-s 

9.85  x 106 

kg-m/ rad 

Ch 

40,500 

kg-m/rad 

Tf 

0.01 

s 

» 


* 
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Azimuth  Axis  Vehicle  Model 

The  azimuth  axis  vehicle  model,  shown  in  Figure  2-3,  is 
the  same  as  the  elevation  model,  except  for  the  values  of  the  moments 
of  Inertia  and  friction.  Angular  positon  is  denoted  by  ip  rather  than  6, 
as  in  the  elevation  model.  A list  of  the  vehicle  azimuth  parameters 
is  contained  in  Table  2-2, 


Table  2-2,  Azimuth  Vehicle  Parameters 


Value 

Units 

8.3  x 104 

kg-m-s 

7.75 

kg-m-s 

3140 

. 2 
kg-m-s 

1.84  x 104 

kg-m-s^ 

9.4  x 106 

kg-m/rad 
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2.2.3  Friction 


Three  types  of  friction,  viscous,  coulomb,  end  st let Ion, 
are  considered  here.  Actuators  generally  have  sigtificant  quantities  of 
all  three  types.  In  addition,  the  turret  and  the  gun  trunnion  may  have 
one  or  more  of  the  friction  types.  The  three  types  of  friction  are 
described  below. 


Viscous  Friction 

Viscous  friction  is  proportional  to  the  speed  of  the  gun 
relative  to  the  hull. 


T - D (0  - 6.  ) 

v v g h 


(2-1) 


Coulomb  Friction 

Coulomb  friction  is  constant  in  magnitude.  The  sign,  or 
direction,  changes  when  the  direction  of  motor  rotation  changes. 


T - D sign  (e  - 8.) 
c c g n 


(2-2) 


Stlctlon  Friction 

Stiction  occurs  each  time  the  motor  starts  or  reverses 
direction.  It  rapidly  decays  to  zero.  The  sign  changes  when  the  direc- 
tion of  rotation  reverses.  Stiction  friction  is  modeled  by  the  follow- 
ing equation: 


Tf  8 • 

3 T- sign  (0  - 0,  ) 

sl+T-8  6 ' g h 


(2-3) 


Total  friction  is  the  sum  of  the  three  types  of  friction  for  both  axes: 


T.  * T + T + T 
f V c s 


(2-4) 


When  compliance  is  added  to  the  model,  the  actuator  and 
trunnion  frictions  have  different  equations.  When  compliance  is  not 
used,  the  equations  are  the  same,  and  the  actuator  and  trunnion  friction 
can  be  combined  as  is  done  here.  The  friction  values  given  in  equa- 
tions (2-1)  through  (2-4)  are  the  combined  friction  of  the  actuator  and 
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trunnion,  for  elevation.  For  the  azimuth  axis,  equation  (2-4)  represents 
the  combined  friction  of  the  actuator  and  the  turret. 

The  block  diagram  used  for  simulating  friction  is  shown  in 

Figure  2-4. 

An  ideal  circuit  for  friction  simulation  is  shown  in  Fig- 
ure 2-5.  In  this  diagram,  the  breakaway  torque  AT  is  the  minimum  actu- 
ator torque  that  will  result  in  actuator  movement.  When  the  magnitude 
of  Ta  is  less  than  AT  and  the  gun  is  at  rest,  the  friction  torqje  exactly 
equals  the  actuator  torque.  At  this  point,  no  torque  is  applied  to  the 
gun.  For  this  region  of  operation,  the  friction  is  called  static  friction. 


Tf  • Ta  (Tf  - static  friction) 


When  the  actuator  torque  exceeds  the  breakaway  torque,  the  circuit 
switches  to  running  friction,  and  Tf  is  given  by  equation  (2-4).  A 
torque  is  then  applied  to  the  gun,  and  gun  motion  starts.  As  long  as 
the  gun  speed  exceeds  the  threshold  value  A6,  the  circuit  stays  in  the 
running  friction  mode  regardless  of  the  value  of  T^.  The  parameters  in 
the  friction  model  must  be  adjusted  so  that  at  the  instant  of  switching 
from  static  to  running  friction,  the  two  friction  values  are  equal. 

2.3  SERVOVALVE  AND  MOTOR  MODELS 

The  servovalve  and  motor  for  both  the  elevation  and  azimuth  systems 
are  modeled  as  shown  in  Figure  2-6.  This  model  includes  a pressure  feed- 
back servovalve  with  nonlinear  flow  dynamics.  The  voltage  applied  to  the 
servovalve  is  the  sum  of  c and  Pf.  The  quantity  e is  the  output  of  the 
control  law  as  described  in  Section  3.  The  quantity  Pf  is  the  output 
of  the  pressure  feedback  network.  The  servovalve  dynamics  Ga  are  repre- 
sented by  a first  order  lag,  i.e.. 


G 


a 


1 

T S + 1 

n 


The  servovalve  gain  is  Ka.  In  the  simulation,  the  deadband  was  assumed 
zero  while  the  valve  area  Ay  was  limited  to  Ay  mav.  The  servovalve 
flow  rate  Qy  is  given  by: 


% 

% 


for  A > 0 
v 


for  A < 0 
v 
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Figure  2-6.  Servovalve  and  Motor  Simulation  Block  Diagram 


where 

C ” flow  coefficient 

P ■ supply  pressure 
s 

p - motor  pressure 

. a 

The  rate  of  change  of  motor  pressure  Pm  is  given  by: 
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where 


B * bulk  Modulus 


V ■ volume 


• displacement  flow  rate,  i.e. 

*m  • '5,  - V 

and  where  RDm  is  the  motor  displacement  reflected  to  the  gun.  For  a 
piston  actuator,  RD^  Is  replaced  by  LAp.  The  Model  is  otherwise  unchanged. 
Leakage  flow  Q*  Is  assuned  proportional  to  the  notor  pressure  and  is  thus 
given  by: 


Qt  - L ?i 


The  resultant  actuator  torque  T^  is  then  given  by: 


T - — - P 
AA  100  ■ 


The  dynamic  pressure  feedback  network  Kp  Gp  is  utilized  for  stabilizing 
and  Inproving  the  servovalve/notor  response.  The  compensator  Gp  is 
used  to  wash  out  the  pressure  feedback  at  low  frequency  or  at  steady 
state.  Pressure  feedback  without  the  washout  would  reduce  the  static 
stiffness.  The  function  Gp  can  either  be  provided  by  a servovalve  de- 
signed to  have  dynamic  pressure  feedback,  or  by  an  electrical  network. 
For  each  axis,  the  gain  Kp  is  a system  parameter  that  was  varied  to 
obtain  the  desired  system  response.  The  pressure  feedback  transfer 
function  is  given  by: 


K G - K t~t~ 

P P P 1 + Tp  » 


Table  2-3  contains  the  parameters  of  the  servovalve  and  motor  for  both 
the  elevation  and  azimuth  systems. 


Table  2-3.  Hardware  Parameter* 


Parameter 


Azimuth 


Linear  Model 


3T/3n 

3T/3e 


Value 

Unite 

10,550 

kg/c*2 

955 

ca/a^/kg/ca2 

Variable 

kg-e 

Variable 

kg-* 

Variable 

kg—a 

Negligible 

kg—2 

0.116 

ca/s-kg/ cm“ 

210  (nominal) 

kg/ca2 

0.07 

a 

1 

ca2/(rad/s) 

0.008 

a 

0.22 

cm2 

2970  (initial  studies) 

3 

cm  /rad 

1311  (optimum  value) 

ca3/rad 

8*5 

c.3 

0.45 

cm2 

1639 

3 1 

a*  /rad 

412 

c3 

2.08  x 104 

kg-nw-s 

6.25  x 104 

kg-m/rad 

0.016 

i 

a 
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2.4  SENSOR  MODEL  DESCRIPTIONS 

2.4.1  Cenaral 

A set  of  five  different  ee nsors  was  selected  by  Cha  Con- 
tract ins  Offiears  for  study.  All  flva  sensors  war*  llkaly  candidates 
for  this  application,  and  It  was  daalred  to  aatablish  tha  parforaanca 
characteristics  of  aach  sansor  on  a ralativa  and  absoluta  baala.  In  an 
approprlata  real la tic  environment.  Proper tlaa  of  tha  electric  rata 
tyros,  tha  General  Electric  hydraulic  rate  sensors,  and  tha  Honeywell 
laminar  vortex  sensor  vara  provided  to  Bendix  by  tha  Contracting  Officers. 
The  modal  of  tha  Airaearch  pneumatic  accelerometer  was  provided  directly 
to  Bendix  by  Airaearch.  The  model  of  tha  Bendix  integrating  acceler- 
ometer wea  obtained  from  Bendix  personnel.  For  aach  of  the  three  rate 
sensors,  a configuration  was  used  in  which  a sensor  was  placed  at  the 
gun  and  in  the  hull.  For  the  two  accelerometers,  a sensor  was  placed 
only  at  the  gun. 

Tha  models  for  tha  sensors  are  presented  in  terms  of 
transfer  functions.  All  sensor  models  have  provisions  for  adding  thres- 
hold and  varying  the  sensor  gains.  All  nominal  sensor  gains  ware  con- 
sidered unity.  Tha  sensor  gain  has  no  effect  on  response  to  coamand 
inputs  or  terrain  disturbances.  Gain  distribution  does,  however,  affect 
the  raaponse  to  noise,  and  drift  due  to  offsets.  Noise  and  drift  were 
not  studied  in  detail  in  this  phase.  When  these  effects  are  studied, 
the  actual  sensor  gains  must  be  used. 

2.4.2  Electric  Rate  Gyro  Model 

The  electric  rate  gyros  which  measure  the  angular  gun 
rate  and  hull  rate  can  be  represented  by  a second  order  trensfer  func- 
tion as  given  by: 


1 + K.+  1 2 

W 2 


4 + 2i  , + 1 .2  8h 

u 2 

n w 

n 


where 


u - 157  rad/s 

XX 

c • 0.7 
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2.4.3  Hydraulic  Rate  Sensor  Modal 

The  hydraulic  raCa  aanaor  ucillxad  in  ehia  acudy  ia  aanu- 
facCurad  by  tha  Canaral  Elactric  Corporation.  Tha  dynamic a of  thla  aanaor 
can  ba  rapraaantad  by  a aacond  ordar  tranafar  function  aa  follows: 


1 ♦ 


it. 


J6« 


1 


$ 


h 


whara 

C - 0.707 
and 

<■>  ■ 40  x 2»  rad/s 
0 

The  following  are  additional  charactariatica  of  thia  aanaor: 

Scale  factor:  0.13  pal/deg/a 

Mas a unbalance  drift:  2 deg/s/g  at  1800  pal 

Noise:  0.6  deg/a  (peak  to  peak) 

Supply  pressure  drift:  0.3  deg/s/percent  change  in 

supply  pressure 

Scale  factor  drift:  0.8  percent /percent  change  in  supply 

pressure 

2.4.4  Hydraulic  Integrating  Accelerometer  Model 

The  hydraulic  integrating  accelerometer  produced  by  Bendlx 
measures  tha  gun  acceleration  0g.  Its  transfer  function  is  given  by: 


9 

** 


Tb 

1 + Tb  * 
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where 


•ad 


• 2*  x 0.1  rad/s 

Th«  sensed  acceleration  la  converted  to  a rata  algnal  0g  by  the  rata 
confutation  network  ahown  below. 


where  Ta  la  act  equal  to  which  la  the  tlae  conatant  of  the  Integrating 
accalaroaeter. 

If  the  aenaor  break  frequency  varies  due  to  teaperature 
changes  or  other  factors,  the  sensed  rate  signal  becoaaa 


Tb  8 

1 + Tb  * 


Tb  l1+\9 

T.  y * Tb  ■ 


e 

s 


Above  the  break  frequency,  the  algnal  la  not  affected  by  the  change  In 
break  frequency.  At  lower  frequencies,  there  Is  a gain  change  propor- 
tional to  the  change  in  break  frequency. 

2.4.5  Laminar  Vortex  Rate  Sensor  Model 

The  transfer  function  for  the  laminar  vortex  sensor  pro- 
duced by  Honeywell  is  a first  order  lag  with  transport  delay,  l.e.. 
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where 


• 0.01  • (transport  daisy) 
t a 0.002  s (time  constant) 

2.4.6  Pnaunatlc  Accelerometer  Model 

Tha  basic  pneumatic  accelerometer  produced  by  Garrett 
Airsearch  has  a second  order  transfer  function,  with  lead-lag  compensa- 
tion, in  a closed  loop  configuration.  The  block  diagram  of  the  sensor 
is  shown  in  Figure  2-7  where 

u - 14.8  rad/s 
n 

c a 0.6 

2 

The  accelerometer  aaturates  at  0.7  rad/s  . 

The  rate  computation  network  which  transforms  the  sensed 
acceleration  into  a rate  is  a simple  integrator. 

2.5  CONTROLLER  MODELS 

The  four  controller  model  variations  developed  for  this  study  are 
described  in  this  section,  and  are  listed  below. 

(1)  Rate  command  with  two  rate  sensors 

(2)  Rate  command  with  single  acceleration  sensor 

(3)  Position  command  with  two  rate  sensors 

(4)  Position  command  with  single  acceleration  sensor 

Two  of  the  variations  considered  are  rate  command  systems  where  the  gunner 
cotmaands  angular  rates.  One  of  the  rate  systems  utilizes  two  rate  sensors 
(gun  and  hull)  while  the  other  utilizes  a single  acceleration  sensor  (gun). 
Two  of  the  models  are  position  control  where  the  gunner  commands  an  angu- 
lar position.  Again,  one  of  these  systems  utilizes  two  rate  sensors 


Figure  2-7.  Block  Diagram  of  the  Pneumatic  Accelerometer 
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while  the  other  utilises  e single  acceleration  eenaor.  Block  dlagraaw 
of  each  of  cheee  concrollera  ere  shown  in  Figures  2-6  through  2-11.  Pro- 
portional end  integral  control  la  utilized  along  with  compensation  networks 
in  the  forward  end  feedback  pethe  end  on  the  input  in  the  poaltlon  command 
system.  The  controller  utilising  e single  acceleration  sensor  requires  e 
rate  computation  network  to  transform  the  sensed  acceleration  into  a rate. 
These  rate  computation  networks  are  defined  with  the  acceleration  sensor 
models  in  Section  2.4.4  end  2.4.6. 
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SECTION  3 

CONTROL  SYSTEMS  APPROACH 


3.1  GENERAL 

General  requirements  of  the  M6QA1  tank  gun  stabilization  system 
ara  as  follows: 

(1)  To  rapidly  attain  a target  by  means  of  gunner  commands  after 
the  target  has  been  sighted  by  the  gunner. 

(2)  To  keep  the  gun  on  target  after  the  gun  has  been  aimed  by  the 
gunner  in  spite  of  hull  motions. 

In  order  to  satisfy  criterion  (1),  it  is  necessary  that  the  response 
of  the  system  be  sufficiently  fast  so  that  the  time  from  sighting  the 
target  until  firing  of  the  first  round  be  minimum.  For  satisfying  cri- 
terion (2),  it  is  necessary  to  minimize  the  effects  of  vehicle  hull  rates 
on  the  tracking  error  of  the  stabilization  system. 

Since  the  design  of  a control  system  for  a given  application  is 
subject  to  all  performance  criteria,  the  criteria  for  the  M60A1  gun 
stabilization  system  are  described  in  Section  3.2.  Note  that  these 
criteria  are  given  in  terms  of  pointing  requirements  as  well  as  general 
performance.  Hardware  acceleration  performance  criteria  are  also  given. 

The  control  systems  approach  utilized  in  this  study  for  achieving 
the  specified  performance  criteria  is  described  in  detail  in  Sections  3.3.1 
and  3.3.2,  for  the  rate  and  the  position  command  inputs,  respectively. 

It  was  specified  at  the  start  of  the  program  that  contemporary  conven- 
tional control  system  techniques  would  be  used  to  achieve  the  desired 
performance.  If  conventional  techniques  could  not  achieve  the  performance 
goals,  more  complex  optimal  procedures  would  be  utilized.  At  this  point 
in  the  development  of  the  complete  stabilization  system,  conventional 
methods  of  synthesis  appear  satisfactory. 

3.2  PERFORMANCE  CRITERIA 

The  following  performance  criteria  were  utilized  in  the  control 
system  design  and  hardware  selection  procedure  for  the  stabilization 
system. 

(1)  Tracking  Error  Requirements 

The  requirement  for  the  tracking  error  is  that  it  be  within 
a circle  of  0.5  mil  diameter  67  percent  of  the  time.  Data  to  verify 
conformance  with  this  requirement  will  not  be  obtained  in  Phase  I. 
Therefore,  a goal  was  set  for  the  tracking  error  to  be  0.5  mil  peak  to 
peak  or  less  for  the  HITPRO  bump  course  at  8 mph. 
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(2)  Frequency  Response  and  Stability  Requirements 

The  frequency  response  goal  was  15  Hz  bandpass.  Bandpass  is 
defined  as  the  frequency  at  which  the  amplitude  ratio  is  >3  dB. 


Cain  margin  and  phase  margin  are  two  criteria  which  are  fre- 


quently  used  to  indicate  the  system  stability.  In  this  study,  the 
following  definitions  are  used  to  establish  stability  goals: 

<•) 

Cain  Margin 

■ 

1.0  - (open  loop  amplitude  ratio  at 
180  deg  open  loop  phase  lag) 

(b) 

Phase  Mai^ ... 

“ 

180  deg  - (open  loop  phase  lag  at 
gain  crossover) 

(c) 

Gain  Crossover 

“ 

frequency  at  which  the  open  loop 
amplitude  ratio  is  1.0  (or  zero  dB) 

A gain  or  phase  margin  of  zero  indicates  an  unstable  system.  The  greater 
the  gain  or  phase  margin,  the  less  likely  it  is  that  instability  will 
occur  due  to  variation  in  control  system  components.  The  goals  are  a 
gain  margin  of  0.5  or  greater  and  a phase  margin  of  35  deg  or  greater. 

(3)  Acceleration  and  Speed  Capability  of  Actuators 

The  goal  for  the  elevation  actuator  system  is  to  accelerate 
to  a speed  of  60  deg/s  in  10  deg  of  gun  rotation,  starting  from  rest. 

The  goal  for  the  azimuth  actuator  system  is  to  accelerate  to  90  deg/s 
in  45  deg  of  turret  rotation,  starting  from  rest. 

(4)  Gain  Variation 

The  control  system  shall  be  designed  to  minimize  the  effects 
of  rate  sensor  gain  variation  and  friction. 

3.3  CONTROL  LAW  DERIVATIONS 

The  control  synthesis  approach  utilized  for  this  study  is  outlined 
in  the  following  two  sections,  for  the  rate  command  input  (rate  control 
system)  and  the  position  command  input  (position  control  system),  respec- 
tively. For  the  rate  control  system,  a hull  rate  cancellation  technique 
and  a lead-lag  stabilization  technique  are  presented.  The  position 
control  system  is  mathematically  compared  to  the  rate  control  system. 

3.3.1  Rate  Control  System 

A simplified  block  diagram  of  a rate  command  stabilization 
system' is  shown  in  Figure  3-1.  The  inputs  are  the  command  rate  9s#and 
the  hull  rate  Gh.  The  output  gun  rate  0g,  which  is  compared  with  9S 
to  generate  an  error,  is  processed  by  the  controller  Gj_.  An  additional 
output  is  the  gun  attitude  angle  6g. 


*,  COMMAND  RATE 
GUN  RATE 


Figure  3-1.  Simplified  Block  Diagram  of  a Rate  Command 
Stabilization  System 


The  objectives  of  the  control  system  are  to  have  the  gun 
rate  follow  the  command  rate  and  to  have  the  gun  angle  9g  stabilized 
with  respect  to  hull  rate  inputs. 

Using  proportional  plus  Integral  control,  it  follows  that: 


G 


1 


+ 


8 


where 


■ proportional  gain 

■ integral  gain 

By  use  of  integral  control,  the  steady  state  tracking  error  TE  - 0S  - 0R 
is  a constant  for  command  rate  and  hull  rate  inputs.  Increasing  the 
gain  Ki  will  decrease  the  tracking  error.  Stability  considerations, 
however,  dictate  a practical  limit  on  K^.  In  order  to  handle  this  prob- 
lem in  an  alternate  way,  two  methods  were  devised  to  reduce  tracking 
errors  due  to  vehicle  motions.  These  two  methods  are  described  in  de- 
tail in  Section  5.'!. 
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The  first  of  these  methods  is  a hull  rate  cancellation 
technique.  The  hull  motion  8jj  is  sensed  and  is  then  used  in  the  control 
laws.  The  control  signal  output  c is  therefore  given  by: 


c 


where 


■ Kj  ■ gain  of  the  hull  motion  effect  on  the  gun  rate 

The  second  method  uses  a single  gun  sensor  and  does  not 
use  a hull  sensor.  The  tracking  error  is  reduced  by  increasing  the 
loop  gain  and  introducing  a lead-lag  network  to  maintain  stability. 
The  control  signal  output  in  this  case  is  given  by: 
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lr  8 
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where 


Position  Control  System 

The  block  diagram  for  position  control  with  proportional 
plus  Integral  control  is  shown  in  Figure  3-2.  The  gunner  generates  an 
angular  positon  command  8S  rather  than  a rate  command  8S.  The  position 
command  signal  is  differentiated  and  compared  with  the  gun  rate  signal 
to  obtain  the  rate  error  signal.  Since  perfect  differentiation  cannot 
be  achieved,  a first  order  lag  wich  time  constant  is  included  in  the 
Input  compensation. 

In  addition,  the  gun  rate  signal  is  integrated  and  compared 
with  the  input  position  command  to  obtain  the  position  error  signal. 
Potentiometer  feedback  was  not  used  as  this  would  have  resulted  in  in- 
creased sensitivity  to  hull  sensor  gain  errors,  and  increased  drift  rate. 

The  position  command  system  can  be  shown  to  be  mathematically 
identical  to  the  rate  system  except  for  the  first  order  lag  T£  in  the 
input  compensation.  (This  is  discussed  in  more  detail  in  Sections  5.2.2 
and  6.3.)  On  this  basis,  the  tracking  error  for  hull  rate  inputs  is 
identical  for  the  rate  and  the  position  control  systems.  In  order  to 


*2'  ‘1 
3.3.2 
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per fora  a meaningful  comparison  of  the  two  control  system  concepts,  a 
rate  command  Input  was  assumed  for  both.  Under  these  conditions,  for 
ramp  position  comaand  Inputs,  the  tracking  error  Is  also  the  same  for 
both  the  rate  and  the  position  systems,  as  long  as  the  input  lag  t*  is 
zero.  If  T£  is  large  enough  to  affect  the  system  response,  or  if 
input  compensation  is  not  used,  response  of  the  position  control  system 
to  rate  inputs  will  be  slower.  However,  one  advantage  with  the  position 
control  system  is  that  a position  command  signal  can  be  used.  This  re- 
sults in  considerably  faster  response  to  the  process  of  aiming  the  gun 
at  a target  for  the  position  control  system,  even  if  input  compensation 
is  not  used. 


Figure  3-2.  Simplified  Block  Diagram  of  a Position  Command 
Stabilization  System 
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SECTION  4 

HARDWARE  SELECTION  PROCEDURE 


4. 1 GENERAL 

For  Che  purpose  of  this  contract,  a study  was  performed  for  select- 
ing the  appropriate  actuator  hardware  in  each  of  the  two  axes  to  meet  the 
following  set  of  specifications  set  forth  by  the  Contracting  Officer's 
Representatives : 

Axis  Slew  Rate  Requirement 

Azimuth  90  deg/s  in  45  deg 

Elevation  60  deg/s  in  10  deg 


These  specifications  were  not  previously  required  for  the  M60A1 
actuators.  A mathematical  model  of  the  hardware  was  thus  derived  (Sec- 
tion 4.3),  and  from  this  model  the  appropriate  servovalve,  accumulator 
volume,  transmission  gear  ratio,  and  motor  displacement  were  selected. 
The  selection  results  for  the  elevation  and  axlmuth  axes  are  given  in 
Sections  4.4  and  4.5,  respectively.  A total  pump  flow  capacity  of 
30  gpm  was  assumed  in  this  study. 

Upon  completion  of  the  study,  the  hardware  was  jointly  selected 
by  Bendix  and  personnel  from  Rock  Island  Arsenal,  for  Inclusion  in  the 
subsequent  computer  study. 

4.2  HARDWARE  CONSIDERATIONS  AND  ASSUMPTIONS 

An  actuator  is  considered  optimal  if  it  achieves  a maximum  turret 
or  gun  acceleration.  If  the  actuator  displacement  is  too  small,  the 
actuator  torque  will  not  be  great  enough  to  rapidly  accelerate  the  load 
inertia.  If  the  actuator  displacement  is  too  large,  the  pump  will  not 
provide  sufficient  flow  to  reach  the  desired  speed.  It  follows  that 
an  optimum  actuator  displacement  exists  which  is  small  enough  to  obtain 
the  desired  speed  with  the  available  pump  flow,  and  large  enough  to  ob- 
tain rapid  acceleration.  Likewise,  the  servovalve  flow  area  should  be 
large  enough  so  that  significant  pressure  drop  does  not  occur  across 
the  servovalve  during  acceleration.  A model  for  the  actuator  system  is 
illustrated  in  Figure  4-1. 

For  the  purpose  of  this  analysis,  it  was  assumed  that  a variable 
delivery  constant  pressure  pump  would  be  used.  The  schematic  for  this 
type  of  pump  is  shown  in  Figure  4-2.  The  pump  servovalve  passes  fluid 
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either  co  or  froa  Che  control  piston,  depending  on  the  sign  of  the  pres- 
sure error.  If  the  pressure  is  higher  then  the  set  pressure,  the  vslve 
pssses  fluid  froa  the  Bain  circuit  into  the  stroking  piston  cylinder, 
thus  decreasing  the  delivery  of  the  pump  so  that  the  pressure  will  be 
reduced.  On  the  other  hand,  if  the  pressure  is  low,  the  vslve  passes 
fluid  to  the  pump  case  and  the  control  spring  forces  the  pump  displace- 
ment to  increase. 

Thus,  as  the  flow  deaand  changes,  the  punp  displacement  automati- 
cally changes  to  maintain  the  supply  pressure  constant,  within  the  flow 
range  of  the  pump. 

The  following  hardware  and  fluid  properties  were  neglected,  since 
they  would  not  affect  the  choice  of  the  optimum  actuator  system. 

• Pump  dynamics 

• Fluid  leakage 

• Fluid  compressibility 

• Friction 

The  pump  dynamics  and  the  compressibility  lag  are  generally  suf- 
ficiently fast  to  have  negligible  effect  on  acceleration.  Use  of  an  ac- 
cumulator slows  the  pump  response.  This  effect  was  neglected  for  this 
study.  The  addition  of  pump  dynamics  would  add  considerably  to  the  com- 
plexity of  the  simulation. 

The  turret  and  gun  acceleration  achieved  in  practice  will  be  a 
little  less  than  the  values  obtained  in  the  study  due  to  the  effects  of 
friction  and  leakage. 

The  parameters  which  must  be  optimized  to  obtain  the  best  perform- 
ance are  the  product  of  transmission  ratio  and  motor  displacement  RDm, 
the  maximum  servovalve  flow  area  (Ay  max) , and  accimilator  volume  Va. 

The  product  M)m  is  the  motor  displacement  reflected  to  the  turret  or  gun, 
respectively. 

4.3  MATHEMATICAL  MODEL  OF  THE  ACTUATOR  SYSTEM 

The  pump  model  assumed  here  provides  that  pump  flow  Qp  equal  to 
that  of  the  motor  Qg  below  the  maximum  pump,  flow  in  accordance  with: 


v*. 


for  < 20  gpa 


% 


20  gpm 


for  >_  20  gpa 


(4-1) 


The  maximum  pump  flow  for  a single  axis  was  taken  to  be  20  gpm.  How- 
ever, a 30  gpm  pump  should  be  used  to  provide  the  required  flow  when 
both  axes  are  activated  simultaneously. 
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It  follows  for  the  axlauth  axis  that  equating  the  aotor  flow  to 
the  servovalve  flow  yields: 


(4-2) 


Squaring  equation  (4-2)  and  solving  for  P^: 


(4-3) 


The  turret  acceleration  is: 


*t"T5oi;  p. 
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The  accuaulator  is  shown  in  Figure  4-1.  The  volume  of  gas  under 
the  diaphragm  is  initially  charged  to  the  supply  pressure.  The  volume 
above  the  diaphragm  is  filled  with  hydraulic  fluid.  When  the  hydraulic 
fluid  is  at  the  desired  supply  pressure,  approximately  one-half  of  the 
accuaulator  volume  is  filled  with  gas.  As  long  as  the  servovalve  flow 
does  not  exceed  the  puap  flow  capability,  the  pump  will  maintain  the 
desired  supply  pressure,  and  the  hydraulic  fluid  and  gas  volumes  will 
not  change.  When  the  valve  flow  exceeds  the  pump  flow,  hydraulic  fluid 
trill  flow  froa  the  accuaulator,  the  gas  volume  will  expand,  and  the 
supply  pressure  will  drop.  This  process  is  adiabatic  and  it  follows 
that: 


V.k  ■ c 


(4-5) 


Differentiating  with  respect  to  time: 

kP  V k_1  V + V k P - 0 

s a a as 
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Solving  for  Pt 


. kP  V 

P - -1! 

rs  V 


(4-6) 


Tha  race  of  change  of  gaa  volume  equals  che  difference  between 
the  motor  flow  and  pump  flow. 


Va  “ ^ “ % • when  > % 


(4-7) 


¥_  » 0 


, when  Q,  < Qp 


For  a hydraulic  motor,  the  dlaplacement  flow  la  given  by: 


*t 


The  gaa  volutae  In  the  accumulator  and  the  supply  pressure  are: 


7 - V ♦ / V dt 

a o J a 


Vq  - steady  state  gas  volume  In  accumulator  at  3000  psl 


P.  * P.o  + k 


Equations  (4-1),  (4-3),  (4-4),  and  (4-7)  through  (4-10)  were  used 
for  the  computer  study  to  determine  the  actuator  system  sizing.  The 
same  set  of  equations  was  used  for  the  elevation  axis,  with  only  a 
change  In  parameters  and  a substitution  of  6g  for 

Table  4-1  lists  the  actuator  system  parameters. 
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Table  4-1.  Parameters  for  Hardware  Selection  Study 


Parameter 

Value 

Units 

C 

100 

ln/s^psl 

J. 

3800 

slugs 

Jc 

22,700 

slugs 

k 

1.4 

- 

Pao 

3000 

psl 

Vo 

one-half 

accumulator 

volume 

in3 

4.4  ELEVATION  AXIS  HARDWARE  SELECTION 

Results  of  the  hardware  selection  procedure  for  the  elevation 
axis  are  summarized  in  Figure  4-3.  Gun  elevation  rate  achieved  at  an 
elevation  angle  displacement  of  10  deg  is  plotted  for  three  different 
servovalve  areas  and  numerous  values  of  the  product  of  transmission 
ratio  R and  motor  displacement  D^.  It  is  evident  from  Figure  4-3  that 
with  a servovalve  area  of  0.034  in2,  the  gun  elevation  rate  exceeds 
60  deg/s  in  10  deg  of  displacement  for  values  of  RDm  between  50  and  110. 
For  maximum  rate,  a value  of  80  for  RD,,,  was  selected  for  elevation. 


Since  servovalves  are  generally  rated  in  terms  of  gpm  flow  at  a 
pressure  drop  of  1000  psl  instead  of  effective  flow  area,  an  appropriate 
conversion  is  achieved  by  use  of  Figure  4-4.  This  figure  also  indicates 
which  commercially  available  valves  can  handle  a given  application.  The 
flow  area  of  0.034  in2  is  thus  equivalent  to  a 20  gpm  servovalve,  and 
the  largest  available  MOOG  series  30  valve  can  be  utilized  for  this 
application.  A 20  gpm  servovalve  was  therefore  selected.  A MOOG 
Series  30  valve  was  also  used  in  an  actuator  system  utilized  in  the  past 
by  the  Chrysler  Corporation  for  the  M60A1  main  gun. 


A time  history  plot  for  the  gun  elevation  angle,  the  elevation 
rate,  and  the  accumulator  gas  volume  is  contained  in  Figure  4-5.  As 
indicated  by  this  figure,  the  accumulator  was  not  used  in  the  first 
10  deg  of  elevation  displacement.  Therefore,  an  accumulator  is  not  re- 
quired for  elevation.  An  elevation  accumulator  may  be  of  value,  how- 
ever, when  both  the  gun  and  turret  are  rotated  simultaneously. 
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9000  PSI.  20  GPM  PUMP  FLOW 

1 GAL.  ACCUMULATOR  USED.  NO  ACCUMULATOR  FLOW 
OCCURRED  DURING  FIRST  10-DEGREE  ROTATION. 

ACCUMULATOR  MAT  HELP  WHEN  AVAILABLE  PUMP 
FLOW  IS  REDUCED  DUE  TO  TURRET  ROTA'.  ION. 

90  GPM  PUMP  REQUIRED  | 

FRICTION  AND  FLOW  LOSSES  NEGLECTED  2 


Figure  4-3.  Hardware  Selection  Results  for  Elevation  Axis 


Figure  4-4.  Servovalve  Conversion  Diagram 
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4.5  AZIMUTH  AXIS  HARDWARE  SELECTION 


Results  of  the  hardware  selection  for  the  azimuth  axis  are  summar- 
ized in  Figure  4-6.  Turret  azimuth  rate  achieved  at  an  azimuth  displace- 
ment angle  of  40  deg  is  plotted  for  three  different  servovalve  areas  and 
numerous  values  of  the  product  RDm.  It  is  evident  from  Figure  4-6  that 
a rate  of  90  deg/s  in  45  deg  displacement  is  not  possible.  It  was  there- 
fore decided  to  choose  the  valve  and  motor  combination  which  maximizes 
rate  at  45  deg  displacement. 

In  accordance  with  Figure  4-6,  a turret  speed  of  70  deg/s  is  obtained 
with  Ay  equal  to  0.07  in?  and  RDm  values  of  100  to  120  in?/rad  using  a 
1 gallon  accumulator.  Therefore,  these  are  the  recommended  values. 

Figure  4-7  demonstrates  the  effect  of  reducing  or  increasing  the 
accumulator  volume  from  1 gallon.  Increasing  the  accumulator  volume  above 

1 gallon  has  no  effect,  while  reducing  the  volume  has  only  a minimal  effect. 

Figure  4-S  shows  the  turret  rotation  required  to  accelerate  to  a 
speed  of  90  deg/s.  If  this  criterion  had  been  used,  an  RDm  of  90  and  a 

2 gal  accumulator  would  have  been  selected.  However,  an  RDm  of  100  ap- 
pears to  be  the  best  overall  choice  when  allowance  is  made  for  torque 
loss  due  to  friction. 

As  indicated  by  Figure  4-4,  the  selected  0.07  in2  servovalve  flow 
area  is  equivalent  to  40  gpm  flow  at  1000  psi  pressure  drop.  This  is 
within  the  range  of  the  MOOG  72  and  73  flow  control  servovalves. 

A time  history  plot  for  the  turret  azimuth  angle,  the  turret  azi- 
muth rate,  and  the  accumulator  gas  volume  is  contained  in  Figure  4-9 
for  the  selected  values  of  RDm  and  Ay. 

A curve  for  each  of  the  azimuth  and  elevation  axes,  which  shows 
the  range  of  values  available  for  R and  Dg,,  is  contained  in  Figure  4-10. 

The  product  RDm  was  determined  to  be  100  and  80  for  the  azimuth  and 
elevation  axes  respectively. 
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SECTION  5 

ANALYTICAL  SYSTEM  ANALYSIS 


5.1  GENERAL 

This  section  contains  an  extensive  analytical  analysis  of  the  rate 
and  positon  control  systems  which  was  performed  prior  to  the  simulation 
analysis.  By  means  of  this  study,  it  was  possible  to  gain  valuable  in- 
sight into  the  behavior  of  the  system  and  to  derive  starting  values  for 
the  rate,  proportional,  and  integral  gains  of  the  control  law.  In  addi- 
tion, two  methods  were  derived  to  reduce  system  tracking  errors.  An 
analytical  sensor  error  study  and  an  outline  of  a coupling  model  are 
also  presented. 

For  the  purpose  of  this  analytical  study,  the  control  system  was 
linearized  and  simplified,  where  possible,  while  retaining  basic  system 
functions.  This  simplification  of  the  system  is  described  in  Sec- 
tions 5.2.1  and  5,2.2  for  the  rate  and  the  position  control  concepts, 
respectively.  Transfer  functions  for  the  system  described  here  showed 
how  the  systems  should  be  designed  to  achieve  the  desired  performance. 

Initial  gain  values  for  the  nonlinear  computer  models  were  also 
determined  from  the  linear  study.  These  gain  values  were  varied,  and 
compensation  was  used  to  achieve  the  desired  performance  with  the  non- 
linear flow  equations,  sensor  dynamics,  and  other  nonlinear  effects  and 
lags. 


5.2  LINEAR  SYSTEM  DERIVATION 

5.2.1  Rate  Control  System 

This  section  contains  the  derivation  of  the  linearized  rate 
control  system  shown  in  Figure  5-1.  The  following  nonlinearities  were 
omitted  in  order  to  simplify  the  transfer  function: 

(1)  Friction 

(2)  Sensor  dynamics 

(3)  Compliance 

Since  pressure  feedback  dynamics  have  negligible  phase 
shift  at  the  natural  frequency,  this  effect  can  be  represented  by  a 
constant  as  shown  in  Figure  5-1. 


5-1 


FOR  READOUT  ONLY 


are  given  by: 


The  nonlinear  flow  and  torque  equations  from  Figure  2-6 


“v  V' ' “•  <8«  ' ^ + * a 


(5-1) 


T.  *RD  P 
A ms 


(5-2) 


Linearizing  by  partial  differentiation  yields: 


f P - P 

C 


CA 


vo 


2 V2  (P8  » P0)  “ 


AP  - RD  s (A9  - A0.)  + LAP 

- g h d 


(5-3) 


Substituting  for  Pm  from  equation  (5-2)  and  rearranging  terms  gives: 


'P  - P 

C _ / -ft-s-  ° AA 


AT.  + 


SAT. 


v 2RD  (P  - P ) A 2BRD  "“‘A 
m s o m 


+ RD  s (A0  - A0.  ) 

m g n 


(5-4) 


where 


Q - CA 
o vo 


J 


P±T “ + 2L  (Ps  - Po> 


Rearranging  terms,  and  defining  new  parameters  gives: 


•rrt^-IK-v] 


(5-5) 
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where 


c 


T 

e 


II 
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3T 
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Assuming 


2CRD  (P,  - p ) 

m ns  o 

CA  + 2L  V2  (P  - P ) 
vo  v s o 

■ 0,  gun  acceleration  is  given  by: 


6 

S 


(5-6) 


The  control  lav  developed  in  Section  3.2  is: 


c 


(5-7) 


Upon  solving  equations  (5-4)  through  (5-7),  the  system  transfer  function 
is: 
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where 
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The  steady  state  T^/Q^  gain  is: 


e 
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3T 
K — 

**1  3c 

3T 

l + KPfr 


(5-12) 


The  response  to  terrain  input  is: 


e 


21  31 

3n  " n 3e 

K — 

i 3c 


1 

G 


(5-13) 


where 

G - cubic  denominator  of  the  transfer  function 
c 

It  follows  that  the  gun  motion,  or  tracking  error,  for  a terrain  input 
is  zero  if 


v . 3T^3n 

Ti  3T/3e 


However,  since  the  torque  speed  slope  is  nonlinear,  perfect  cancellation 
over  the  complete  range  of  operation  may  not  be  possible. 

5.2.2  Position  Control  System 

This  section  contains  a derivation  of  the  linearized  ele- 
vation position  control  system  shown  in  Figure  5-2.  As  indicated  in  the 
figure,  the  inertial  gun  position  feedback  signal  is  obtained  by  integrating 
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th«  gun  rat*  sensor  signal.  A potentiometer  was  not  used  for  tha  position 
faadbaek  because  It  would  have  bean  necessary  to  add  th*  hull  position  to 
tha  potentiometer  signal  to  obtain  the  inertial  position  signal.  This 
would  have  wade  the  system  very  sensitive  to  hull  sensor  gain  variation. 

The  control  law  for  thla  system  is  given  by: 


. • <e.  - e,>  ♦ kc  (5,  5„ 


Sine*  a position  coamand  input  will  generally  be  used  for  th*  position 
command  system,  the  coamand  rate  signal  6g  will  not  be  readily  available. 
This  signal  can  be  omitted,  or  it  can  be  obtained  by  pseudo  differenti- 
ation of  the  input  eomaand  signal  as  shown  below. 


e_  - 


1 + S 


Including  the  term  results  in  faster  response  to  command  inputs. 

Other  equations  for  the  position  control  concept  take  the 
same  form  as  those  for  the  rate  control  concept. 

It  follows  therefore  that  the  position  control  transfer 
function  is  given  by: 


s 


1 ♦ 


1 + T, 


1 + 
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2 3 

+ 25_  + _±_ 

*4  *4 


The  steady  state  response  to  a hull  rate  step  input  la: 
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The  transfer  fuse  cions  for  hull  rata  sad  disturbance  torque  inputs  are 
the  ease  aa  for  the  rate  control  spates.  Therefore,  the  response  to 
these  Inputs  is  identical  for  the  two  systi 


The  response  to  identical  rate  cn— nd  inputs  will  be  the 
ease  for  both  systems  if  the  pseudo  differentiator  lag  t*  equals  zero. 
This  condition  was  uasd  in  the  simulation  for  shoving  the  equivalence 
of  the  two  systems,  as  discussed  in  Section  6.3.  Since  position  command 
inputs  will  generally  be  used  with  the  position  control  system,  the 
response  time  to  rotate  the  gun  to  a commanded  angle  will  be  shorter 
whether  or  not  the  input  rate  algnal  0,  is  used. 

3*2.3  Linear  Model  Parameters 

Por  model  development  and  preliminary  studies,  the  actuator 
displacement  and  volume  under  compression  were  obtained  from  the  bibliog- 
raphy. The  parameters  developed  here  corresponded  to  the  elevation  axis: 


«D 


LA  - 38.4  in.  x 4.72  in  x 16.37 


in' 


where 

RD^  - 2970  cm3/rad 
V - 895  cm3 

The  supply  pressure  and  maximum  servovalve  area  were  taken  to  be 

P#  - 210  kg/cm2 


Ay  - 0.2  cm 


The  linearized  actuator  parameters  are: 


il . (RD  ,2  Vila 

3n  100  '“V  Q 
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V 

B 


II 

ac 


The  effective  servovalve  area  Ay0  and  pressure  P0  are 
typically  taken  to  be  equal  to  one-half  the  maximum  servovalve  area  and 
one-half  supply  pressure  for  linear  analysis.  Thus,  the  linearized 
servovalve  parameters  are: 

/P.  " ?o  /210  - 105 

Qo  <*vo  V 2 955  x °*1  V 2 

giving 

Q - 695  cm3/s 
o 


Also 


fjf  - jfio  (2970) 2 HI  - 2.66  x 104  kg-o-s 


895 


m - 


c 10,550  * 695 


0.013  s 


||  - 0.02  x 2970  x — - 6.25  x 104  kg-m/rad 


Generally,  the  linearized  parameters  must  be  varied  some- 
what to  obtain  the  best  correlation  with  the  nonlinear  system.  To 
achieve  this  correlation,  the  following  values  were  used: 

H ■ 2.08  x 104  kg-m/s 

x ■ 0.016  s 
c 

■ 6.25  x 104  kg-m/rad 

0 £ 
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The  Initial  gain  values  for  K*,  Kf  and  Kp  ware  determined 
as  follows.  The  bandpass  goal  was  IS  Hz.  Thus,  a bandpass  of  16  Hz 
was  selected  for  the  linear  system,  to  allow  some  performance  margin. 
From  equation  (5-9) 


v . 1 / ,2  3T  | 

Kr  3T/3e  |Tc  JuNS  " 3n  j 


and  for 


* 2*  x 16  • 100  rad/s 


it  follows  that: 


t - r (0.016  x 527  x 1002 *  - 2.08  x 104 *)  - 1.0 

6.25  x 10* 


From  aquation  (5-10) 


, _ L . 3T/3n 
Ci  %S  [Kr  3T7fr 


In  this  expression,  B should  equal  0.35  for  the  best  linear  system  re- 
sponse, and  a should  equal  0.7.  Thus: 


2 OR 

K±  - 0.35  x 100  1.0  + - 47 


Equation  (5-11)  specifies  that: 


Kp  " 3T737  (aTc  “NS  - X) 


Upon  substituting  numerical  values,  it  follows  that: 


K - 7 (0.016  x 0.7  x 100  - 1)  - 2 x 10"6 

P 6.25  x 104 


u 


f 


5.2*4  Linear  System  Response 

Baaed  on  the  transfer  functions  presented  in  Sections  5.2,1 
and  5.2.2,  a digital  computer  solution  for  time  and  frequency  response 
vas  obtained.  Using  ideal  sensors  and  the  parameters  for  the  elevation 
axis,  the  resulting  time  response  is  illustrated  in  Figure  5-3.  Mote 
that  the  time  response  exhibits  a 7 percent  overshoot.  A solution  for 
system  frequency  response  is  illustrated  in  Figure  5-4.  This  figure 
shows  that  a bandpass  of  18  Hz  was  achieved. 


5.3  TRACKING  ERROR  REDUCTION 

This  section  deals  with  the  method  of  reducing  the  gun  tracking 
error  due  to  motions  of  the  hull.  Referring  to  Figure  5-5,  which  is  a 
simplified  block  diagram  of  the  rate  command  stabilization  system  with 
proportional  and  integral, control,  assume  a constant, hull  rate  input 
and  no  command  input  (6S  ■ 0).  In  steady  state,  9„  » 0 and  Qv  - Qq. 
Then: 


e 


3T/3n  • 
VtlTc  °h 


Since  the  tracking  error  TE  is  given  by: 


TE  - 8 - 8 

• g 


e - Kr  <eg  - eg)  + Kt  <e8  - eg> 


then  the  steady  state  tracking  error  is: 


TE  - 


e 


3T/3n  • 

K±  3T/3e  h 


It  is  desired  to  reduce  this  tracking  error  so  that  the  gun  is  stabi- 
lized with  respect  to  hull  motions. 

Two  methods  can  be  used  to  reduce  the  tracking  error:  The  first 

method  is  a cancellation  technique.  Equation  (5-13)  shows  that  hull 
motion  will  cause  no  gun  motion,  and  thus  will  result  in  no  tracking 
error,  if 


„ 3T/3n 
ti  3T/3e 


! 


i 

! 


{ 

\ 
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e CONTROLLER  OUTPUT 
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0"  MOTOR  PLOW 
ta  ACTUATOR  TORQUE 


figura  5-5,  Simplified  Block  Diagram  of  Che  Rate  Comnand  System 


Therefore,  hull  motion  can  be  effectively  cancelled  by  using  the  output 
signal  of  an  inertial  hull  sensor  with  gain  in  the  control  law.  Com- 
plete cancellation  of  hull  motion  coupling  for  the  actual  system,  how- 
ever, cannot  be  obtained  due  to  nonlinearities  and  compressibility.  In 
addition,  gain  variations  due  to  such  factors  as  aging  and  temperature 
will  cause  some  tracking  error.  The  cancellation  technique,  however, 
is  vary  effective  in  reducing  gun  rot ion. 

The  second  method  uses  only  a .single  gun  sensor.  No  vehicle  motion 
signal  input  is  provided  for  this  system.  The  tracking  error  for  the 
system  will  be  inversely  proportional  to  the  error  signal  amplification 
and,  therefore,  a very  high  gain  must  be  used.  A lag-lead  network  is 
required  to  lower  the  high  frequency  gain  for  stability.  The  design 
technique  is  to  first  adjust  the  gains  for  the  stable  operation  with 
the  desired  response.  The  gain  required  to  lower  the  tracking  error  to 
an  acceptable  level  is  then  determined.  The  lag-lead  network  is  designed 
to  maintain  the  high  frequency  gains  at  the  original  values  for  stable 
operation,  while  a higher  gain  is  obtained  at  low  frequency  for  low  track- 
ing error. 

The  higher  gain  and  lag-lead  netvork  can  also  be  used  with  the  can- 
cellation technique  to  reduce  the  sensitivity  to  variation  of  some  of 
the  gains. 
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5.4  SENSOR  ERROR  ANALYSIS 


5.4.1  General 


The  effect  of  sensor  errors  on  the  control  and  stabiliza- 
tion systea  is  the  subject  of  this  section.  Sensor  errors  considered 
here  ate  gain  errors  and  offset  errors.  Only  the  rate  comnand  system 
is  analyzed  in  this  section.  It  is  assumed  that  in  both  the  rate  and 
position  comnand  systems,  the  gunner  will  zero  the  tracking  error  with 
or  without  sensor  errors.  On  that  basis,  sensor  errors  have  an  effect 
only  on  stabilizing  the  gun  in  the  presence  of  hull  motions.  As  described 
in  Section  5.2.2,  the  response  to  hull  motions  is  identical  for  both  the 
rate  and  the  position  command  systems. 


utilizing: 


The  analysis  presented  here  will  consider  a system 


(1)  Two  rate  sensors  (Section  5.4.2) 


(2)  A single  acceleration  sensor  in  the  gun  axis 
(Section  5.4.3) 


5.4.2  Two  Rate  Sensors 


Figure  5-6  is  a simplified  block  diagram  of  the  rate  con- 
trol system  with  two  rate  sensors.  The  transfer  function  of  the  sensors 
Is  assumed  to  be  Y(s)  with  unity  gain.  For  zero  sensor  gain  error, 

Kbh  “ Kbg  " 1*  The  sensor  offset  errors  are  and  tg. 

A hull  sensor  gain  error  has  the  same  effect  as  a change 
in  Kj,.  The  effects  of  changes  in  Kj,  were  investigated  using  the  simula- 
tion and  are  described  in  Section  6. 


A hull  sensor  offset  error  will  give  a constant  steady 
gtatea tracking  error.  This  is  determined  as  follows:  Assuming  that 


Bjj  ■ 0 0g  ■ 0 (in  steady  state),  then  e ■ 0 and  -K*  eg  » Kh  ej,. 


Since  the  tracking  error  TE  ■ -6g, 


TC  - rt  eh 


Typical  gain  values  are  fy,  ■ 0.15  and  * 50.  Therefore,  the  tracking 
error  in  radians  for  an  offset  error  in  rad/s  is  TE  * 0.003  ej,. 


A gun  sensor  gain  error  will  have  little  effect  on  the 
stabilization  system  since  it  only  changes  the  loop  gain  and  is  like 
changing  Kr  and  K<.  The  effect  of#a  gun  sensor  offset  error  Eg  is  analyzed 
as  follows:  Assuming^g  « 0,  and  0j,  - constant,  then  in  steady  state. 


’gs 


TE 


0.  Therefore,  0g  ■ ~eg*  Since  the  tracking  error  is  given  by 


e 


g,  it  follows  that  TE  - ft g dt.  The  tracking  error  for  a gun 
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Figure  5-6.  Simplified  Block  Diagram  of  the  Rate  Control  System 
with  Two  Rate  Sensors 


sensor  offset  will  therefore  increase  with  time.  As  a result,  it  is 
mandatory  to  minimize  the  gun  sensor  offset  error  or  devise  a method 
of  eliminating  its  effect.  Section  5.5  discusses  a method  of  eliminating 
sensor  offset  errors. 

5.4.3  Single  Acceleration  Sensor 

A simplified  block  diagram  of  a rate  control  system  with 
a single  acceleration  sensor  in  the  gun  axis  is  presented  in  Figure  5-7. 
Both  an  integrating  and  a pneumatic  accelerometer  are  shown,  along  with 
the  corresponding  rate  computation  networks. 

For  both  of  these  acceleration  sensors,  a gain  error  (i.e., 
Kbg  5*  1)  1135  the  same  effect  as  the  rate  sensors.  The  loop  gain  changes 
with  the  gain  error,  but  the  influence  on  the  tracking  error  is  small. 
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Figure  5-7.  Simplified  Block  Diagram  of  the  Rate  Control  System 
srith  a Single  Acceleration  Sensor 

For  the  Integrating  accelerometer,  the  effect  of  an  offset 
error  is  synthesized  as  follows:  In  steady  state*  0gs  - 0.  Therefore, 


• 1 b “ 

9 _ - 0 - 1 + — — 2— r 0 + e 

gs  i s t.  s + 1 g g 


It  follows  that 


6 ■ -s  e - 

g g 


g Tfe  S g 


Since  Che  cracking  error  is  TE  - -0 


8* 


TE  - 


i‘> ic  * 


-*•  dt  dc 
Tb 


The  cracking  error  will  Cherefore  increase  vich  cime  for  a consCanC  off- 
aeC  error,  ice  effecC,  however,  can  be  minimized  by  increasing  Che  Cime 
conscanc  tj,.  For  Che  pneumadc  acceleromeCer,  Che  effect  of  an  offseC 
error  is  decermined  in  a similar  manner.  In  aCeady  scace,  0g8  - 0 so 
ChaC 


1 * 

0 ■ — (0 

8 gS 


The  Cracking  error  is  cherefore 

IE  - V JJ‘* dt 

Again  Che  cracking  error  will  increase  wich  Cime  for  a constant  sensor 
offset  error. 

5.5  METHOD  OF  ELIMINATING  THE  EFFECTS  OF 

INTEGRATOR  DRIFT  AND  SENSOR  OFFSETS 

The  sensor  error  analysis  in  Section  5.4  shows  that  gun  sensor 
offsets  cause  Che  gun  to  drift,  and  hence  Che  cracking  error  Co  increase 
with  time.  A method  has  been  developed  which  automatically  nulls  out 
Che  effects  of  integrator  drift  and  sensor  offsets.  This  method  is 
illustrated  in  Figure  5-8.  The  sensor  and  integrator  shown  are  parts 
of  the  original  system.  The  track/store  amplifier  operates  such  that 
the  output  tracks  (equals)  the  input  when  in  the  Crack  mode,  and  stores 
(holds)  the  output  when  in  the  store  mode. 

The  amplifier  is  ordinarily  in  the  store  mode  so  that  Che  output 
of  the  sensor  8„s  (gun  rate  in  this  case)  is  integrated  to  give  the  gun 
angle.  At  the  time  it  is  desired  to  null  the  sensor  offset  eg  and  the 
integrator  drift  e<j,  the  track/store  amplifier  is  put  into  the  track 
mode. # This  must  be  accomplished  during  a period  when  the  actual  gun 
rate  0g  is  zero.  Under  these  conditions,  when  the  network  reaches  steady 
state,  its  output  is  zero,  effectively  nulling  the  offset  and  drift. 
Switching  the  amplifier  to  the  store  mode  returns  the  system  to  normal 
operation  except  that  now  the  drift  and  offset  are  nulled  out. 
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Figure  5-8.  Method  of  Eliminating  the  Effects  of  Integrator  Drift 
end  Sensor  Offsets 


Two  of  these  networks  are  required  for  an  acceleration  sensor 
system  while  only  one  is  required  for  a rate  sensor  system.  In  addi- 
tion, the  network  is  required  only  for  the  gun  sensor  since  an  offset 
error  in  the  hull  sensor  results  in  a small  tracking  error  which  does 
not  increase  with  time. 

In  order  to  implement  this  type  of  drift  and  offset  elimination, 
the  gunner  can  be  provided  with  a pushbutton  switch  which  activates 
ell  the  nulling  networks.  A timer  could  also  be  utilized  to  automati- 
cally return  the  system  to  normal  operation  after  the  nulling  process 
is  completed. 

5.6  EFFECT  OF  GYROS  ON  RESPONSE 

A block  diagram  for  analysis  of  the  effect  of  rate  sensor  dynamics 
on  response  is  presented  in  Figure  5-9.  The  transfer  function  for  a 
command  input  is 
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where 


0^  ■ denominator  of  gyro  transfer  function 

K#  " gain  of  actuator  and  gun  transfer  function 

■ denominator  of  the  actuator  and  gun  transfer  function 

The  gyro  dynamics,  represented  by  Dg,  add  a lead  term  to  the 
transfer  function,  which  trill  increase  the  overshoot  to  a step  input 
command.  The  gyro  dynamics  also  add  phase  lag  to  the  open  loop  sys- 
tem, and  thus  will  reduce  the  damping.  Therefore,  compensation  must 
be  used  to  reduce  the  overshoot  and  phase  lag  resulting  from  the  sensor 
dynamics. 

5.7  FORWARD  VERSUS  FEEDBACK  COMPENSATION 

Lead-lag  compensation  must  be  used  to  reduce  the  overshoot  and 
phase  lag  resulting  from  the  sensor  dynamics.  In  this  section,  a com- 
parison is  made  between  using  this  type  of  compensation  in  the  forward 
or  in  the  feedback  path. 

The  block  diagram  in  Figure  5-10  shows  a system  open  loop  transfer 
function  with  a gain  term  Ks  in  the  numerator,  dynamic  terms  Ds  in  the 
denominator,  and  with  feedback  compensation  Nc/Dc. 

The  dosed  loop  transfer  function  is 
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The  block  diagram  in  Figure  5-11  shows  the  same  system  with  for- 
ward compensation.  The  closed  loop  response  is 
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The  response  with  feedback  compensation  has  a lead  term  equal  to 
the  compensation  denominator  Dc.  The  response  with  forward  compensation 
has  a lead  term  equal  to  the  compensation  numerator  Nc.  Generally,  the 
smaller  the  lead  time  constant,  the  less  the  overshoot  will  be.  There- 
fore, when  lead-lag  compensation  is  used,  it  should  be  in  the  feedback 
path.  When  lag-lead  compensation  is  used,  it  should  be  in  the  forward 
path. 
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5.8  SYSTEM  STIFFNESS 

Stiffness  is  the  reciprocal  of  the  response  to  a disturbance 
torque.  The  static  stiffness  can  be  considered  as  the  equivalent  spring 
rate  of  the  system.  The  tracking  error  resulting  from  a constant  dis- 
turbance torque  equals  the  disturbance  torque  divided  by  the  static 
stiffness.  The  stiffness  transfer  function  is 


T.  1 + a.  s + a,  s2  + a,  s3  + a.  s4 

Ji  m k 1 1 2 4 

9g  ° (1  + Tp  .)  (1  + tc  .)  + Kp  Tp  ff  s2 
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The  static  stiffness  is  given  by 
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The  stiffness  should  be  high  to  minimize  the  effects  of  disturbance 
torques.  Disturbance  torques  can  result  from  gun  bending  mode  vibrations, 
coupling,  and  bumps  in  the  road.  The  static  stiffness  can  be  increased 
by  increasing  the  system  integral  gain  K^.  When  is  Increased,  a low 
frequency  dipole  must  be  used  in  the  forward  path  to  lower  the  gain  at 
high  frequency  as  required  for  stability. 

The  analysis  in  Section  5.9  shows  the  effect  of  stiffness  on  track- 
ing errors  resulting  from  coupling. 
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5.9  COUPLING 

The  azimuth  and  elevation  atabllizatlon  aysteaa  are  coupled  by  both 
gyroscopic  and  hull  roll  rate  effects.  Thla  aectlon  dlacusaea  thla 
coupling  and  methods  of  minimizing  their  effect. 

One  cause  of  coupling  between  the  azimuth  and  elevation  axes  is 
gyroscopic  moments.  According  to  the  lav  for  a gyroscope,  if  a mass  is 
rotating  about  the  x axis  with  angular  momentum  Hx,  simultaneous  rotation 
about  the  y axis  will  result  in  a torque  about  the  z axis.  Thus, 


Ts  - • Hx  (5-14) 

where 

H - J w (5-15) 

x xx  ' ' 


Applying  this  law  to  the  elevation  axla,  with  the  gun  pointing  straight 
ahead,  gives: 


gc 


gy  yt  yt 


Jgx  *h 


(5-16) 


where 

J • inertia  of  the  gun  about  the  roll  axis 
gx 

“ inertia  of  the  gun  about  the  yaw  axis 

Thus,  when  turret  rotation  and  roll  occur  simultaneously,  the  resulting 
torque  Tgc  is  applied  to  the  gun. 

The  gun  moments  of  inertia  about  the  yaw  and  roll  axes  are  not 
known.  A rough  estimate  of  the  inertia  about  the  yaw  axis  was  calcu- 
lated by  multiplying  the  gun  mass  by  the  square  of  the  estimated  distance 
from  the  gun's  center  of  gravity  to  the  center  of  turret  rotation. 


J - 1.52  x 100  - 225  kg-m-s2 

gy 

The  gun  moment  of  inertia  about  the  roll  axis  is  assumed  to  be  negligible. 
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The  peak  roll  rate,  from  the  HITPRO  program,  la  0.24  rad/s.  The 
peak  turret  rate  la  assumed  to  be  1.571  rad/s. 

The  torque  applied  to  the  gun  by  the  gyroscopic  effect  is  thus 
Tgc  " ♦h  • Jgy  * " °*24  * 225  x \ - 85  kg-« 


The  tracking  error  resulting  from  this  torque  is  Tgc  divided  by 
the  stiffness.  The  stiffness  curve  determined  from  the  nonlinear  com- 
puter simulation  of  the  rate  control  system  with  electric  gyro  is  shown 
in  Figure  5-12.  The  stiffness  at  0.6  Hz,  the  fundamental  frequency  of 
the  bump  course,  is  5.6  z 10^  kg-m/rad.  Thus 


TE  - — =•  - 0.00015  rad  (0.15  mil) 

5.6  x 10J 


When  the  gun  is  not  pointing  straight  ahead,  the  coupling  equa- 
tion becomes 


Tgc  " (Jgx  * Jgy>  (*h  c°8  *t  + ®h  ata  V *t  (5“17) 


" (Hj> 


Figure  5-12.  Dynamic  Stiffness  Curve 
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Similarly,  for  the  azimuth  axis,  a torque  Ttc  is  applied  to  the 
turret  when  the  vehicle  rolls  while  the  gun  is  being  rotated  about  the 
pitch  axis,  and  when  vehicle  pitch  and  roll  occur  simultaneously. 


where 


Jtz  6h  + 6h  * Jtx  *ta  + *h  ' Jg  8g 


Jtx  " inertia  of  turret  about  the  roll  axis 

J • inertia  of  turret  about  the  pitch  axis 
cz 

The  turret  inertias  are  assumed  to  be  equal  about  all  three  axes. 
The  torque  applied  to  the  turret,  using  hull  rates  from  HITPRO,  is 


T • 2 J l,  0.  + ♦,  J 0 - 2 x 3140  x 0.24  x 0.18  + 0.24  x 527  x 1.05 
tc  t h h h g g 

* 408  kg-ra 


The  dynamic  stiffness  curve  was  not  run  for  azimuth.  However, 
azimuth  stiffness  should  be  about  six  times  higher  than  for  elevation, 
due  to  the  higher  inertia  of  the  turret.  The  higher  stiffness  in  azi- 
muth was  indicated  by  the  smaller  tracking  errors  on  the  bump  course 
with  this  axis.  The  estimated  tracking  error  resulting  from  gyro  effect 
is  thus 


6 x 5.6  x 10J 


0.00012  rad  (0.12  mil) 


The  estimated  tracking  errors  are  within  the  goal  of  0.25  mil,  but 
are  large  enough  to  be  significant.  The  calculated  values  could  be  con- 
siderably in  error  if  the  estimated  inertias  are  significantly  in  error. 

Two  methods  can  be  used  to  reduce  the  tracking  errors  resulting 
from  coupling.  The  first  method  is  by  increasing  the  stiffness.  Stiff- 
ness is  proportional  to  the  integral  gain  K^,  as  shown  in  Section  5.8. 
Therefore,  increasing  Ki  will  reduce  the  coupling  error.  It  is  also 
important  to  tune  the  system  parameters  and  compensation  so  that  the 
stiffness  remains  high  over  the  frequency  band  of  coupling. 
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The  second  method  uses  a hull  roll  rate  sensor  to  cancel  the  coup- 
ling effect.  The  coupling  error  can  probably  be  reduced  to  an  acceptable 
level  by  increasing  the  stiffness.  Therefore,  the  second  method  was  not 
studied  in  detail. 

A second  coupling  effect,  for  elevation,  is  the  result  of  hull 
pitch  angle,  as  seen  by  the  gun,  varying  with  turret  angle.  When  the 
turret  angle  is  zero,  the  torque  acting  on  the  gun  is  given  by 


V 


where  Dg^  is  the  trunnion  friction  and  3T/3n  is  the  slope  of  the  actu- 
ator torque-speed  curve.  The  significance  of  the  3T/3n  term  will  now 
be  discussed.  Figure  5-13  shows  the  gun  with  a piston  actuator,  and 
with  the  servovalve  represented  by  orifice  area  Ay.  If  there  were  no 
oil  in  the  cylinder,  and  no  trunnion  friction,  hull  pitch  motion  would 
not  apply  any  torque  to  the  gun.  The  gun  inertia  would  then  maintain 
the  gun  orientation  fixed  in  space.  However,  when  there  is  oil  in  the 
cylinder,  any  motion  of  the  hull  relative  to  the  gun  results  in  oil 
being  pumped  through  orifice  Ay.  This  flow  causes  a differential  pres- 
sure to  act  on  the  piston,  and  thus  applies  a torque  to  the  gun.  The 
ratio  of  the  torque  to  the  relative  velocity  is  the  term  3T/3n. 


Figure  5-13 


Schematic  for  Coupling  of  Hull  Pitch  Rate  to  Gun 


The  equation  for  friction  coupling  becomes  considerably  mor»  com- 
plex when  vehicle  roll  is  considered. 

A totque  also  results  from  hull  acceleration.  This  torque  is 
proportional  to  the  gun  mass  unbalance  Lg  Mg  and  the  distance  between 
the  hull  centers  of  rotation  and  the  trunnion. 


T,2  • L,  ",  «p  9h  “*  + *«■  V 


The  mass  unbalance  of  the  gun  is  small.  Therefore,  this  term  will  prob- 
ably be  negligible. 

A second  coupling  for  azimuth  results  from  the  turret  mass  unbalance 
with  respect  to  the  turret  center  of  rotation  L[  Mt.  When  the  turret 
angle  is  zero,  hull  roll  motion  applys  a torque  to  the  turret,  but  hull 
pitch  has  no  effect.  When  the  turret  angle  is  90  deg,  hull  pitch  motion 
applies  a torque  to  the  turret,  and  hull  roll  has  no  effect.  The  equa- 
tion for  the  torque  applied  to  the  turret  by  these  effects  is 


Tt  ■ (sin  9h  sin  ^ + sin  ^ cos  * ) 


5.10  SUMMARY  OF  ANALYTICAL  SYSTEMS  STUDY 

Following  is  a summary  of  the  results  obtained  and  methods  of 
synthesis  accomplished  by  means  of  the  analytical  systems  study.  Some 
of  these  items  proved  to  be  ir  valuable  in  the  subsequent  computer  analy- 
sis of  the  system. 

(1)  A linear  model  of  the  stabilization  system  was  derived  which 
aided  in  the  design  of  the  control  system  and  allowed  the 
establishment  of  an  initial  set  of  control  gains. 

(2)  For  terrain  inputs,  the  responses  of  the  rate  and  position 
command  systems  are  identical. 

(3)  Two  methods  can  be  used  to  reduce  the  tracking  error  result- 
ing from  hull  motion.  These  are: 

(a)  Cancellation  of  hull  rate  using  hull  rate  sensor 

(b)  Lag-lead  compensation  with  high  loop  gain 

(4)  A sensor  error  analysis  restricted  to  gain  and  offset  errors 
showed  that  only  gun  sensor  offset  errors  are  significant 
and  cause  the  tracking  error  to  increase  with  time. 

(5)  A method  was  developed  for  eliminating  the  effects  of  sensor 
offset  errors  and  integrator  drift. 
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(6)  Sensor  dynamics  Increase  overshoot  and  reduce  stability. 

(7)  Lead-lag  compensation  will  be  required  if  the  gun  rate  sensor 
phase  lag  is  significant. 

(8)  When  a lead-lag  configuration  is  used,  it  should  be  in  the 
feedback  path. 

(9)  The  integral  control  gain  K*  should  be  high  to  maximize 
system  stiffness  and  minimize  the  effects  of  disturbance 
torques. 

(10)  When  vehicle  roll  rates  occur,  significant  coupling  between 
the  azimuth  and  elevation  axes  takes  place  in  the  form  of 
gyroscopic  moments  and  other  effects. 
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SECTION  6 

SIMULATION  ANALYSIS  RESULTS 


6.1  GENERAL 

This  section  contains  a detailed  description  of  the  analog  and 
digital  computer  results  obtained  in  this  study  for  the  gun  stabiliza- 
tion system.  The  objectives  of  the  computer  study  are  outlined  here. 

The  organization  of  the  relevant  computer  program  as  well  as  a summary 
of  the  types  of  inputs  and  outputs  are  discussed  in  Section  6.2.  Sec- 
tion 6.3  compares  the  rate  and  position  command  systems.  The  results 
of  « preliminary  computer  analysis  evaluating  the  effects  of  nonlinear 
valve  flow,  vehicle  hull  dynamics,  and  hydraulic  fluid  compressibility 
are  described  in  Section  6.4.  The  effects  of  these  nonlinearities 
needed  to  be  evaluated  in  order  to  judge  whether  or  not  they  would  be 
included  in  the  detailed  simulation  models.  A detailed  evaluation  of  the 
complete  stabilization  system  and  the  five  candidate  sensors  are  contained 
in  Sections  6.5  and  6.6. 

The  following  is  a summary  of  the  objectives  considered  in  this  sy- 
tems  study. 

(1)  Evaluate  the  effects  of  the  nonlinear  valve  flow  equations 
versus  linearized  flow,  and  the  effect  of  hull  dynamics. 

(2)  Derive  stabilization  system  configurations  for  each  of  the 
sensors,  and  define  gains  and  compensation  networks  required 
to  meet  the  performance  goals  with  each  system. 

(3)  Investigate  the  effects  of  sensor  gain  variation,  sensor  dead- 
band, and  actuator  friction  on  tracking  error. 

(4)  Perform  detailed  evaluation  of  the  five  prospective  sensors 
for  this  application. 

(5)  Determine  the  system  time  response  for  step  and  sinusoidal 
inputs. 

(6)  Determine  the  tracking  error  for  sinusoidal  null  rates  and 
for  bump  course  using  the  HITPRO  digital  computer  program. 

(7)  Determine  the  gain  and  phase  margins  for  each  configuration. 

The  following  sensors  were  analyzed  in  detail. 

(1)  Elevation  Axis 

Two  electric  gyros. 

Two  hydraulic  rate  sensors  (GE) . 

Two  laminar  vortex  rate  sensors  (Honeywell). 

One  integrating  accelerometer  (Bendix). 
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(2)  Azimuth  Axis 

Two  electric  gyros 

One  pneumatic  accelerometer  (AlResearch) 

Rate  sensors  can  be  used  In  either  one  (gun)  or  two  (gun  and  hull) 
sensor  configurations.  Only  the  rate  sensor  configuration  utilizing  two 
sensors  was  analyzed  in  this  study. 

Accelerometers  can  be  used  in  only  the  one  sensor  configuration. 

Since  the  accelerometer  signal  is  integrated  to  obtain  the  rate  signal, 
an  accelerometer  for  sensing  hull  rate  would  have  an  open  loop  integrator 
which  would  result  in  excessive  drift  rates.  On  the  other  hand,  the  gun 
sensor  integrator  is  within  a closed  control  loop.  Under  those  conditions, 
integrator  drift  will  have  less  effect. 

6.2  SIMULATION  APPROACH 

The  organization  of  the  computer  programming  approach  utilized  for 
the  analog,  hybrid,  and  digital  computer  programs  is  described  in  this 
section.  The  types  of  inputs,  organization  of  the  analog  computer  boards, 
analog  computer  switch  logic  utilized,  and  the  hybrid  computer  inputs 
are  also  described.  Detailed  listings  of  all  digital  computer  programs 
are  contained  in  Appendix  A.  Analog  and  hybrid  program  wiring  diagrams 
are  presented  in  Appendix  B. 

Digital  computer  solutions  were  obtained  for  each  system  configura- 
tion as  checks  on  the  analog  computer  simulations.  The  following  inputs 
were  generally  used  for  the  digital  computer  simulations. 

(1)  Command  rate  step  inputs 

(2)  Hull  rate  step  inputs 

(3)  Hull  rate  sinusoidal  inputs 

The  analog  program  of  the  control  system  was  initiated  by  means  of 
a simplified  linear  model.  The  nonlinear  flow  equations,  sensor  dynamics, 
nonlinear  friction,  and  deadband  were  added  one  at  a time,  so  that  their 
effects  could  be  assessed  individually. 

In  order  to  facilitate  the  selection  of  the  various  components  for 
a particular  computer  simulation  run,  automated  computer  switches  were 
used  for  the  following  purposes. 

(1)  Selecting  one  of  the  five  sensors. 

(2)  Selecting  either  a two  sensor  configuration  using  a hull  sensor 
and  a gun  sensor  or  a one  sensor  configuration  using  only  a 
gun  sensor. 

(3)  Selecting  either  forward  or  feedback  path  compenc ation. 
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The  following  Inputs  were  used  for  the  hybrid  computer  study  for 
determining  time  and  frequency  responses. 

(1)  Step  command  rate  input. 

(2)  Sinusoidal  rate  command  input. 

(3)  Step  hull  rate  input. 

(4)  Sinusoidal  hull  rate  input. 

(5)  Simulated  bump  course  from  program  HITPRO. 

The  step  and  sinusoidal  command  inputs  used  were  sufficiently  small 
to  avoid  saturation  of  any  component  of  the  simulation.  A hull  rate  in- 
put in  the  azimuth  axis  as  used  here  is  representative  of  a pivot  steer 
maneuver.  The  digital  program  HITPRO  was  used  to  provide  realistic  hull 
motion  for  the  simulated  bump  course.  To  implement  these  rates  in  the 
simulation,  the  hull  pitch  motion  was  stored  in  the  digital  part  of  the 
hybrid  computer  from  a magnetic  tape  and  was  then  used  for  providing  the 
elevation  hull  rate  inputs.  The  same  rate  input  was  used  for  the  azimuth 
axis  because  no  comparable  azimuth  rates  were  available  from  the  HITPRO 
program.  This  is  a valid  procedure  since  the  power  spectral  densities 
of  the  hull  pitch  and  yaw  motions  are  approximately  the  same. 

The  sinusoidal  hull  rate  input  used  was  of  approximately  the  same 
amplitude  as  the  maximum  hull  rate  observed  for  the  bump  course,  and  of 
approximately  the  same  frequency  as  the  bump  course  fundamental  frequency. 

6.3  COMPARISON  STUDY  OF  RATE  AND  POSITION 

CONTROL  FOR  STABILIZATION 

As  discussed  in  Section  5.2.2,  the  transfer  functions  for  hull  rate 
inputs  are  the  same  for  both  the  rate  and  position  command  systems.  Also, 
the  transfer  functions  for  command  inputs  are  the  same  if,  for  the  posi- 
tion control  system,  the  input  compensation  term  s/(ii  > + 1)  is  a per- 
fect derivative. 

In  order  to  verify  the  above  results,  a computer  study  was  made  to 
compare  responses  using  the  rate  and  position  command  simulation  models. 
Each  of  the  rate  and  position  command  systems  was  modelled  with  equal 
gains  and  compensation  networks  and  with  two  electric  gyro  sensors. 

A perfect  derivative  of  the  position  command  0a  was  used  in  the  position 
system.  The  frequency  response  curves  of  the  two  systems  are  shown  in 
Figures  6-1  and  6-2.  Figure  6-1  verifies  that  the  two  systems  are,  in 
fact,  equivalent. 

In  an  actual  system,  a perfect  derivative  cannot  be  generated; 
instead,  an  approximation  to  it  must  be  generated  by  the  input  compen- 
sation network  of  the  type  s/(x£  s + 1).  Figure  6-2  shows  the  system 
frequency  response  using  a perfect  derivative  and  with  * 0.01  s for 
the  input  compensation  network. 

As  a result  of  this  study,  the  performance  evaluations  of  each  of 
the  five  sensor  configurations  were  conducted  using  the  rate  conmand 
system  only. 
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6.4  PRELIMINARY  SIMULATION  STUDIES 


6.4.1  Effect  of  Nonlinear  Valve  Flow 

It  was  found  by  means  of  the  computer  simulation  that  the 
rate  command  system  (with  ideal  sensors)  exhibited  basically  the  same 
step  response  with  the  linear  flow  models  as  with  the  nonlinear.  When 
the  electric  gyro  dynamics  were  added,  the  nonlinear  model  became  un- 
stable, while  the  linear  model  remained  stable.  This  indicates  that 
the  nonlinear  flow  model  has  less  phase  margin  than  the  linear  model. 

The  nonlinear  flow  model  was  therefore  used  for  all  subsequent  system 
studies. 

6.4.2  Effect  of  Hull  Dynamics 

Hull  dynamics  modelled  in  this  simulation  represent  the 
effect  observed  when  the  actuator  applies  torque  to  the  gun.  At  that 
point  an  equal  and  opposite  torque  is  applied  to  the  hull  and  the  result- 
ing hull  motion  affects  the  gun.  The  simulation  was  used  to  determine 
system  response  to  command  and  hull  rate  inputs,  with  and  without  hull 
dynamics.  It  was  found  that  hull  dynamics  had  a negligible  effect  on 
system  response  and  on  the  tracking  error.  As  a result,  hull  dynamics 
were  not  included  in  the  model  for  subsequent  studies. 

6.4.3  Effects  of  Changes  in  Fluid  Compressibility 

and  Bulk  Modulus 

The  volume  under  compression  V modelled  in  this  simulation 
consists  of  the  high  pressure  fluid  in  the  servovalve,  actuator,  and 
connecting  lines.  The  numerical  value  for  V used  here  was  obtained 
from  the  bibliography.  This  value  seems  large  and  indicates  that  in  the 
past,  long  connecting  lines  were  used.  To  determine  whether  this  signi- 
ficantly affected  system  performance,  a step  response  simulation  was 
performed  with  a 10  percent  reduction  in  the  value  of  V.  A negligible 
change  in  the  step  response  of  the  system  was  observed. 

The  effect  of  a drastic  change  in  the  bulk  modulus  was 
also  determined.  This  effect  represents  the  process  of  air  entering 
the  hydraulic  fluid.  It  was  found  that  reducing  the  fluid  bulk  modulus 
B by  50  percent,  has  a negligible  effect  on  the  system  step  response  and 
stability.  These  results  indicate  that  the  amount  of  pressure  feedback 
used  here  was  great  enough  to  compensate  for  the  large  volume  under 
compression. 

Compressibility  may  have  different  effects  when  compliance 
and  gun  bending  modes  are  added  to  the  model.  These  effects  should 
therefore  be  evaluated  again  in  future  studies  when  bending  modes  and 
compliance  have  been  modelled. 

6.5  EVALUATION  OF  SENSORS  IN  ELEVATION  AXIS 

Sections  6.5  and  6.6  contain  the  evaluation  of  each  of  the  five 
sensors  studied.  Time  and  frequency  response  plots  are  presented  to 
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illustrate  the  results.  The  results  for  the  sensors  In  Section  6.5.1 
through  6.5.4  are  based  on  elevation  (gun)  axis  parameters  while  the 
results  for  the  pneumatic  accelerometer  in  Section  6.6  are  based  on  azi- 
muth axis  parameters.  All  simulation  results  presented  here  are  based 
on  the  rate  command  control  concept  as  pointed  out  in  Section  6.3. 

6.5.1  Evaluation  of  Electric  Rate  Gyros 

A block  diagram  of  the  rate  command  system  with  two  elec- 
tric rate  gyros  is  presented  in  Figure  5-6.  A hull  rate  sensor  is  used 
for  cancellation  of  hull  rate  coupling,  and  a gun  rate  sensor  is  used 
for  closing  the  control  loop.  Lead-lag  compensation  is  used  to  compen- 
sate for  sensor  phase  lag,  and  thus  to  Improve  stability.  Compensation 
is  in  the  feedback  path.  This  is  the  best  location  for  the  compensation 
in  order  to  reduce  the  overshoot  to  command  step  Inputs.  (See  Sec- 
tion 5.7.) 

Using  the  hybrid  computer  simulation,  the  gains  were  ad- 
justed to  obtain  the  desired  frequency  response,  stability,  and  tracking 
error.  The  resulting  gain  values  and  compensation  are: 


K±  - 50  l/s 


Kp  « 2.5  x 10"4  rad/(8-kg-cm2) 

- 0.15 

Feedback  compensat  ion  fr^nnVt8 

1 + 0.0016  s 

The  response  to  a command  step  input  0S  - rad/s  was  ob- 
tained, and  is  shown  in  Figure  6-3.  The  gun  rate  reaches  the  commanded 
speed,  0.05  rad/s,  in  about  0.04  s,  overshoots  about  15  percent,  then 
settles  to  the  commanded  speed  with  no  oscillations.  The  absence  of 
oscillations  indicates  good  damping.  In  the  absence  of  a specific  per- 
formance requirement  for  step  response,  the  speed  of  response  as  well  as 
damping  appear  satisfactory. 
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Figure  6-3 • Response  to  a Command  Step  Input  - Elevation  Rate 
Control  with  Electric  Gyros 


The  error  race  signal  ca  given  by  e8  - e8  - 0ga  la  the 
input  for  Che  ooen  loop  response.  The  open  loop  response  Is  given  by: 


e 

s 


The  variable  6g8  (sensed  gun  rate)  i|  Che  output  signal 
for  the  open  loop  response.  The  variables  6g  and  0g8  are  used  mainly 
for  determining  frequency  response,  gain  margin,  ana  phase  margin. 
These  variables  are  also  used  as  an  aid  in  tuning  the  system. 

The  tracking  error  TE  and  the  sensed  tracking  error  £7 
are  defined  as: 


te  - e - 0 

* g 

e_  * f c dt  ■ 9 - 0 

T J s s gs 

The  variable  TE  originates  at  an  open  loop  integration  and  therefore 
tends  to  drift.  The  variable  ex  is  obtained  from  within  the  closed 
loop  and  generally  does  not  drift.  These  two  signals  are  nearly  iden- 
tical except  when  sensor  deadband  or  a high  frequency  input  are  used. 
The  actuator  torque  T^  is  also  shown. 

The  response  to  a sinusoidal  hull  pitch  rate  is  shown  in 
Figure  6-4.  The  pitch  rate  input  is: 


0.  ■ 0.18  sin  (?ir  x 0.6  t)  rad/s 
n 


This  input  is  approximately  equal  to  the  most  severe  pitch 
rate  experienced  on  the  simulated  bump  course.  As  shown  in  Figure  6-4, 
the  resulting  tracking  error  has  a peak-to-peak  amplitude  of  0.18  mil 
and  thus  meets  the  pointing  performance  specification  of  the  system. 


s • 

The  closed  loop  frequency  response  curves  0g/6g  are  shown 
in  Figure  6-5.  The  open  loop  frequency  response  curves  were  used  for 
determining  the  gain  and  phase  margins.  An  input  signal  amplitude  of 
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0.03  rad/s  peak  to  peak  was  used  for  all  frequency  response  tests.  The 
frequency  response  results  are  as  follows. 


System  bandpass: 

Open  loop  phase  lag 
at  gain  crossover: 

Phase  margin: 

Gain  margin: 


28  Hz 
108  deg 

180  deg  - 108  deg 
0.5 


72  deg 


The  system  gain  margin  was  determined  by  increasing  the  gain  until  in- 
stability resulted.  It  was  found  that  the  system  was  stable  with  the 
gain  doubled,  and  unstable  with  higher  gain.  This  is  equivalent  to  a 
gain  margin  of  0.5. 


The  frequency  response  results  meet  or  exceed  the  per- 
formance criteria  specified  in  Section  3.3. 


Time  response  plots  using  inputs  consisting  of  vehicle 
rates  from  the  simulated  bump  course  of  the  HITPRO  program  are  given 
in  Figures  6-6  through  6-8.  Graphs  summarizing  these  results  are  pre- 
sented in  Figures  6-9  and  6-10.  The  maximum  tracking  error  for  nominal 
conditions  (Kj,  - 0.15)  was  0.18  mil  peak  to  peak. 

Figure  6-9  shows  that  a hull  sensor  gain  error  of  50  per- 
cent increased  the  tracking  error  to  0.4  mil. 


Figure  6-10  shows  that  combined  coulomb  and  stiction  fric- 
tion of  50  kg-m  each,  increased  the  tracking  error  to  0.4  mil. 

The  0.18  mil  tracking  error  observed  with  nominal  condi- 
tions is  well  within  the  0.5  mil  criteria  listed  in  Section  3.2.  Note 
however  that  the  combined  effects  of  sensor  gain  error,  sensor  deadband, 
friction,  noise,  and  cross-coupling  of  axes  will  increase  the  tracking 
error.  The  tracking  error  can  then  be  reduced  by  increasing  the  system 
low  frequency  gain.  This  can  be  accomplished  by  adding  lag-lead  compen- 
sation in  the  forward  path,  of  the  form: 


12 

T1 


1 + T, 


1 + T. 


S 


S 


A graph  demonstrating  dynamic  stiffness  is  presented  in 
Figure  6-11.  These  results  were  obtained  by  applying  a sinusoidal  dis- 
turbance torque,  and  recording  tracking  error  as  the  output  signal. 
Stiffness  is  the  ratio  of  the  torque  amplitude  to  the  tracking  error. 
The  minimum  stiffness  was  found  to  be  111  dB.  The  dynamic  stiffness 
should  be  high  to  minimize  the  effects  of  coupling  and  bending  mode 
vibration. 
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6.5.2  Evaluation  of  Hydraulic  Rate  Sensor 

Figure  5-6  is  a block  diagram  of  the  stabilization  system 
with  the  hydraulic  rate  sensor.  The  basic  configuration  is  the  same  as 
with  two  electric  gyros.  Control  gains  and  feedback  compensation  were 
identical  to  those  used  with  the  electric  gyros. 

Usin';  {he  hybrid  computer,  the  response  of  the  system  to 
a step  rate  command  6S  was  obtained.  As  shown  in  Figure  6-12.  this 
step  response  exhibited  a 20  percent  overshoot.  The  response  time  is 
about  the  same  as  for  the  system  with  two  rate  gyros. 

The  response  to  sinusoidal  hull  motion  of  an  amplitude 
equivalent  to  the  bump  course  is  shown  in  Figure  6-13.  The  tracking 
error  observed  was  0.12  mil  peak  to  peak. 

The  frequency  response  results  shown  in  Figure  6-14  can 
be  summarized  as  follows. 

Bandpass:  13  to  37  Hz 

Gain  margin:  0.53 

* Phase  margin:  50  deg 

Note  that  the  amplitude  ratio  is  down  3 dB  at  13  Hz  and  drops  to  -5  dB 
at  21  Hz.  Modification  of  the  feedback  compensation  network  will  allow 
improvement  of  the  system  bandpass. 

The  response  to  the  HITPRO  simulated  bump  course  is  shown 
in  Figures  6-15  through  6-17.  Figure  6-18  shows  the  maximum  observed 
tracking  error  versus  hull  sensor  gain  error.  A gain  error  of  33  per- 
cent increased  the  maximum  tracking  error  ro  0.38  mil.  The  maximum 
tracking  error  versus  sensor  threshold  or  deadband  is  shown  in  Fig- 
ure 6-19.  A deadband  of  6 mils  increases  the  tracking  error  to  2 mils. 
The  effects  of  combined  deadband  and  friction  effects  are  shown  in  Fig- 
ure 6-20.  A small  amount  of  friction  har  whe  effect  of  reducing  the 
cracking  error.  Larger  amounts  of  fric  ion  increase  the  tracking  error. 
The  combination  of  deadband  and  friction  increases  the  tendency  for 
limit  cycling. 

6.5.3  Evaluation  of  Integrating  Accelerometer 

The  block  diagram  of  the  stabilization  system  with  an  inte- 
grating accelerometer  is  presented  in  Figure  5-7.  Only  a single  sensor, 
mounted  on  the  gun,  is  used.  The  processed  sensor  signal  is  proportional 
to  acceleration  at  frequencies  below  0.1  Hz  and  to  rate  above  0.1  Hz. 

In  the  simulation  model, the  compensation  network  was  designed  to  inte- 
grate the  sensor  signal  below  0.1  Hz  and  thus  to  provide  a rate  signal 
at  all  frequencies. 

Since  there  is  no  hull  sensor,  the  system  gain  must  be 
increased  to  reduce  the  tracking  error.  To  achieve  this,  compensation 
was  used  in  the  forward  path.  The  type  of  compensation  increases  the 
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Figure  6-20.  Effects  of  Combined  Hull  Sensor  Deadband  and  Friction  - 
Hydraulic  Rate  Sensor 

gain  at  low  frequencies  and  decreases  the  gain  at  high  frequencies  as 
Is  required  for  stability.  The  forward  compensation  and  system  gains 
established  using  the  hybrid  computer  simulation  are  as  follows: 

Forward  compensation:  10  * + 0*55^* 

Control  gains:  K*  ■ 250 

Kr  - 7.5 

Kp  ■ 2.5  x 10"4 

The  step  response  of  this  sensor  configuration  to  a command 
rate  Input  Is  shc*n  in  Figure  6-21.  The  response  has  a 60  percent  over- 
shoot. This  overshoot  can  be  reduced  significantly  using  a lead-lag 
feedback  compensation. 

The  response  to  a sinusoidal  hull  input,  equivalent  in 
amplitude  to  the  most  severe  motion  for  the  bump  course  Is  shown  In 
Figure  6-22.  The  tracking  error  observed  was  0.11  mil  peak  to  peak. 

The  response  to  the  actual  HITPR0  bump  course  is  shown  in  Figure  6-23. 
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Friction  - Elevation  Rati 
Control  with  Integrating 


The  maximum  tracking  error  observed  is  0.28  mil  peak  to  peak.  The 
tracking  error  for  the  HITPR0  input  could  be  reduced  by  modifying  the 
compensation. 

Coulomb  and  stiction  friction  values  of  100  kg-m  each 
increased  the  tracking  error  to  0.38  mil,  as  shown  in  Figure  6-24. 

Figure  6-25  shows  the  maximum  tracking  error  for  various  values  of 
friction. 

The  frequency  response  curves  are  shown  in  Figure  6-26. 

The  frequency  response  results  can  be  summarised  as  follows. 

Bandpass:  16  Hz 

Gain  margin:  0.53 

Phase  margin:  45  deg 

6.5.4  Evaluation  of  Laminar  Vortex  Sensor 

A block  diagram  for  the  rate  command  system  with  two  lami- 
nar vortex  sensors  is  shown  in  Figure  5-6.  The  sensor  models  consist 
of  a single  lag  and  a time  delay.  The  time  delay  was  simulated  in  the 
digital  part  of  the  hybrid  computer.  To  achieve  this  for  the  gun  rate 
sensor,  the  gun  rate  signal  from  the  analog  computer  was  fed  into  the 
digital  computer,  delayed  in  time,  and  then  fed  back  to  the  analog  com- 
puter. For  the  hull  sensor,  the  hull  rate  from  the  HITPRO  program  was 
delayed  in  time  in  the  sensor  path,  but  used  without  time  delay  in  the 
remaining  portion  of  the  simulation. 

The  step  response  to  a rate  command  input  is  shown  in  Fig- 
ure 6-27.  A 30  percent  overshoot  is  exhibited. 

The  response  to  the  simulated  bump  course  using  the  HITPRO 
program  is  shown  in  Figures  6-28  and  6-29  for  hull  sensor  gains  of  0.15 
and  0.20,  respectively.  A plot  of  tracking  error  versus  hull  sensor 
gain  Kj,  is  shown  in  Figure  6-30.  From  this  graph,  a gain  of  0.15  was 
selected  as  the  nominal  value  with  the  minimum  tracking  error.  Increas- 
ing Kj,  to  0.20  increased  the  tracking  error  to  0.6  mil.  Decreasing  Kh 
to  0.10  increased  the  tracking  error  to  0.28  mil. 

The  sinusoidal  hull  input  response  was  not  obtained  for 

this  sensor. 

The  frequency  response  curves  are  shown  in  Figure  6-31. 

The  frequency  response  results  are  summarized  as  follows. 

Bandpass:  25  Hz 

Gain  margin:  0.50 

Phase  margin:  50  deg 
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6.6  EVALUATION  OF  SENSORS  IN  AZIMUTH  AXIS 


6.6.1  Evaluation  of  Electric  Rate  Gyros 

The  block  diagram  of  the  azimuth  stabilization  system  with 
two  electric  rate  gyros,  shown  in  Figure  5-6,  is  identical  to  the  block 
diagram  for  the  elevation  axis.  Only  some  of  the  parameters  were  changed 
to  allow  for  hardware  differences  between  exes.  (See  Tables  2-1  and  2-2.) 

Using  the  hybrid  computer  simulation,  the  step  response  of 
the  system  to  a rate  command  input  was  obtained  and  is  shown  in  Fig- 
ure 6-32.  The  observed  overshoot  is  36  percent. 

The  response  to  sinusoidal  hull  motion  equivalent  in  ampli- 
tude to  the  bump  course  is  shown  in  Figure  6-33.  The  observed  tracking 
error  was  0.03  mil  peak-to-peak. 

The  response  of  the  system  using  the  simulated  HITPRO  bump 
course  is  shown  in  Figure  6-34.  The  maximum  tracking  error  is  else 
0.03  mil  peak  to  peak.  The  response  of  the  system  to  the  bump  course 
with  a 20  percent  hull  sensor  gain  error  was  also  investigated  and  the 
results  are  as  shown  in  Figure  6-35.  This  amount  of  gain  error  increased 
the  tracking  error  to  about  0.07  mil. 

The  effects  of  hull  sensor  gain  variations  for  both  the 
bump  course  and  sinusoidal  hull  motions  are  summarized  in  Figure  6-36. 

It  is  evident  that  the  observed  tracking  error  is  approximately  the  same 
for  either  input. 

Additional  simulation  studies  were  conducted  to  determine 
the  response  to  the  bump  course  with  2 and  6 mil/s  sensor  deadband. 

These  results  are  shown  in  Figures  6-37  and  6-38,  respectively. 

The  effects  of  deadband  and  friction  on  the  tracking  error 
are  summarized  in  Figure  6-39.  Coulomb  friction  of  100  kg-m  increased 
the  tracking  error  to  0.5  mil.  The  addition  of  an  equal  amount  of  stic- 
tion  friction  did  not  change  the  tracking  error.  Sensor  deadband  of 
2 mils/s  increased  the  tracking  error  to  0.07  mil.  For  larger  deadband, 
a limit  cycle  occurred,  and  the  tracking  error  was  greatly  increased. 

The  frequency  response  curves  are  shown  in  Figure  5-40. 

The  frequency  response  results  are  numerated  as  follows. 

Bandpass:  21  Hz 

Gain  margin:  0.50 

Phase  margin:  55  deg 

The  response  of  the  system  to  a pivot  steer  maneuver  is 
shown  in  Figure  6-41.  In  a pivot  steer  maneuver,  the  tank  is  turning 
at  its  maximum  rate,  i.e.,  180  deg  in  8 s.  In  practice,  it  ic  desired 
that  the  gun  remain  pointing  at  the  target  during  this  maneuver.  As 
shown  in  Figure  6-41,  there  is  a momentary  tracking  error  of  0.09  mil. 

The  tracking  error  then  decays  to  zero  after  about  0.5  s. 
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Figure  6-34,  Response  to 
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Electric  Gyros 


□ HULL  MOTION  FROM  HITPRO 
DIGITAL  COMPUTER  PROGRAM 

O HULL  MOTION  - 8.18  SIN  2*  x 0.6 1 


HULL  SENSOR  GAIN  VARIATION,  PERCENT  * 

Figure  6-36.  Tracking  Error  Versus  Hull  Sensor  Gain  Error  - Azimuth 
Rate  Control  with  Electric  Gyros 


6.6.2  Evaluation  of  Pneumatic  Accelerometer 

A block  diagram  of  the  stabilization  system  for  the  azimuth 
axis  using  a pneumatic  accelerometer  is  presented  in  Figure  5-7.  The 
accelerometer  is  used  for  sensing  turret  acceleration,  and  the  sensor 
output  signal  is  integrated  to  obtain  the  turret  rate  signal. 

A high  gain  loop  was  used  to  reduce  the  tracking  error. 
Compensation  was  used  in  the  forward  path  for  reducing  the  gain  at  high 
frequency  to  achieve  stability.  In  addition,  feedback  compensation  was 
used  to  reduce  the  effects  of  sensor  phase  lag. 

The  gains  and  compensation  networks  used  in  the  simulation 
were  as  follows. 


K - 10 
r 

K - 2 x 10-4 
P 

Forward  compensation:  10 


Feedback  compensation: 


1 ± 0.05  s 
1 + 0.5  s 

1 + 0.010  s 
1 + 0.002  s 


urea* 


Figure  6-41.  Response  to  Pivot  Steer  Maneuver  - Azimuth  Rata  Control 
with  Electric  Gyros 


The  atep  reaponae  obtained  with  the  hybrid  computer  simulation  la  shown 
In  Figure  6-42.  An  overshoot  of  44  percent  was  observed. 

The  response  to  sinusoidal  hull  wot ion  equivalent  In  ampli- 
tude to  the  maximum  hull  notion  for  the  bump  course  la  shown  In  Fig- 
ure 6-43.  The  observed  tracking  error  was  0.2.  all  peak  to  peak. 

The  response  to  the  HITPRO  bump  course  Is  shown  In  Fig- 
ure 6-44.  The  effect  of  adding  100  kg-m  coulomb  friction  and  100  kg -a 
stictlon  friction  is  shown  In  Figure  6-45.  The  resulting  maximum  tracking 
error  was  again  0.25  mil.  The  friction  caused  an  Increase  In  the  maximum 
tracking  error  for  larger  values  of  friction,  as  shown  In  Figure  6-46. 

The  frequency  response  curves  are  shown  in  Figure  6-47. 

The  frequency  response  results  can  be  summarized  as  follows. 

Bandpass:  15  Hz 

Gain  margin:  0.55 

Phase  margin:  36  deg 
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Figure  6-46.  Maximum  Tracking  Error  Versus  Friction  - Pneumatic 
Acc llerometer 
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6.7  SUMMARY  OF  SIMULATION  RESULTS 

Simulation  of  the  stabilization  system  has  produced  the  following 
results. 

(1)  The  linear  servovalve-actuator  model  provided  a means  of 
estimating  required  system  gains  but  was  not  sufficiently 
accurate  to  he  used  in  the  sensor  evaluation  study.  The 
nonlinear  model  was  therefore  used  f^r  sensor  evaluations. 

(2)  Hull  dynamics  had  a negligible  effect  on  system  response  and 
tracking  error. 

(3)  A summary  of  the  results  of  the  sensor  evaluation  study  is 
presented  in  Table  6-1. 

Table  6-1 (A)  summarizes  the  results  for  the  rate  sensors. 

For  each  of  these  systems,  a rate  sensor  was  mounted  in  the 
gun  and  the  hull  axis. 

Table  6-1 (B)  summarizes  the  results  for  the  accelerometers 
studied.  For  both  of  these  systems,  an  accelerometer  was 
mounted  in  the  gun  axis  only. 

It  is  to  be  noted  that  the  systems  represented  in  the  table 
were  not  fully  optimized  in  this  study.  An  absolute  compari- 
son of  the  sensors  is  therefore  not  entirely  justified. 

As  indicated  by  Table  6-1,  the  performance  criteria  were  met 
for  all  of  the  sensor  configurations. 

(4)  A hull  sensor  gain  error  of  25  percent  can  be  tolerated 
without  exceeding  the  tracking  error  requirement. 

(5)  Coulomb  friction  of  100  kg-m  increases  the  tracking  error  by 
approximately  0.5  mil  peak  to  peak. 
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SECTION  7 
CONCLUSIONS 


A mathematical  model  of  a suitable  stabilization  system  for  the 
M60A1  tank  main  gun  was  formulated  and  programed.  The  model  was  de- 
fined so  as  to  Include  most  of  the  significant  nonlinearities  such  as 
nonlinear  valve  flow  and  hull  dynamics  due  to  gun  motion.  A hybrid  com- 
puter analysis  was  performed  to  determine  the  operating  characteristics 
of  the  stabilization  system,  to  evaluace  prospective  sensors  for  sensing 
gun  and/or  hull  rate,  and  to  determine  whether  a rate  or  position  com- 
mand control  concept  is  preferable  with  respect  to  specified  performance 
criteria. 

An  analytical  study  revealed  that  both  the  rate  and  position  con- 
trol concepts  required  a proportional  plus  integral  control  law  in  order 
to  minimize  the  gun  tracking  error.  It  was  also  shown  that  the  rate  and 
position  concepts  are  equivalent  in  terms  of  nulling  out  the  effects  of 
hull  motions  and  thus  of  stabilizing  the  gun  after  the  target  Is  in  the 
sight.  A computer  analysis  which  followed  verified  this  equivalence. 

In  addition,  it  was  possible  to  show  that  the  effect  of  hull  motions  on 
the  system  can  be  minimized  by  either  a high  control  loop  gain  along 
with  a lead-lag  compensation  network  or  by  using  a hull  sensor  signal 
in  the  control  law. 

The  extensive  computer  simulation  analysis  revealed  several  con- 
clusions in  the  areas  of  stabilization  control  philosophy  and  sensor 
applicability.  In  the  process  of  arriving  at  a full  computer  model  of 
the  system  for  sensor  evaluation,  it  was  found  that  the  effect  of  hull 
dynamics  on  the  gun  was  negligible.  Nonlinear  valve  flow,  however,  was 
found  to  have  a significant  Influence  on  system  performance.  A linear- 
ized flow  model  was  not  sufficiently  accurate  for  use  in  this  study. 

It  was  found  that  all  five  of  the  sensorB  studied  in  this  program 
meet  the  performance  criteria  set  forth  by  the  Contracting  Officer's 
Representative.  In  addition,  this  study  indicates  that  these  criteria 
can  be  mat  by  using  only  a gun  sensor.  If  verified  by  further  studies, 
the  need  for  a corresponding  hull  sensor  may  be  eliminated. 

The  detailed  sensor  st«dy  lamed lately  revealed  that  automatic  off- 
set and  Integrator  drift  nulling  circuits  are  required  when  using  an 
acceleration  sensor.  A method  which  can  be  used  for  this  purpose  is 
described  in  Section  5.5. 

In  addition,  increasing  the  gain  of  the  acceleration  sensor  will 
decrease  the  sensor  offset  effects,  and  hence  the  drift  rate.  More 
generally.  It  was  found  in  the  sensor  stu^y  that  sensor  gain  errors  have 
a small  effect  on  the  tracking  error.  Also,  a combination  of  sensor 
deadband  and  gun  or  turret  friction  will  cause  the  system  to  limit  cycle. 

In  order  to  compensate  for  sensor  phase  lag,  feedback  compensation 
is  required.  Forward  path  compensation  is  desirable  for  obtaining  sta- 
bility with  higher  loop  gains  for  this  system. 
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SECTION  8 
RECOMMENDATIONS 


It  is  recommended  that  in  future  system  studies  and  sensor  evalua- 
tions, additional  system  characteristics  including  realistic  bending 
modes  and  hardware  compliance  model  be  added.  The  effects  of  providing 
an  extensive  model  for  friction  which  separates  the  effects  of  stiction 
coulomb,  and  running  friction  through  switching  logic  should  also  be 
investigated.  Sensor  models  should  be  utilized  allowing  for  sensor  errors 
due  to  noise,  deadband,  gain  variation,  and  offset.  It  is  also  of  im- 
portance on  a complete  system  model  to  include  integrator  drift  and  an 
automatic  nulling  circuit.  When  these  effects  have  been  included  in  the 
model,  a statistical  analysis  of  the  pointing  error  output  data  will  be 
required  in  order  to  afford  more  detailed  and  objective  performance  com- 
parisons. Inclusion  of  the  additional  nonlinearities  mentioned  will  also 
make  it  possible  to  better  optimize  the  system  gains  and  compensation 
networks. 

It  is  also  recommended  that  a gunner  model  be  included  in  the  track- 
ing loop,  in  order  to  arrive  at  additional  data  for  comparing  the  rate 
and  position  command  control  concept.  A display  of  tracking  error  can 
be  developed  using  an  oscilloscope  for  which  a human  operator  can  issue 
realistic  commands  to  the  system  to  simulate  target  tracking.  Inclusion 
of  an  operator  model  in  the  simulation  will  allow  a determination  of  the 
tracking  and  stabilization  capabilities  of  the  system. 

In  addition,  it  is  recommended  that  steps  be  taken  to  verify  the 
results  of  this  simulation  study  by  means  of  field  test  data  with  an 
M60A1  tank.  Additional  insight  into  the  system  operating  characteris- 
tics and  the  effects  of  significant  nonlinearities  could  be  gained.  A 
verification  study  would  consist  of  obtaining  recordings  of  measurements 
of  hull  and  gun  rates  and  achieved  pointing  accuracies  on  various  types 
•»f  terrain,  and  comparing  these  with  the  results  of  this  study. 

An  optimal  control  theory  approach  to  the  stabilization  of  the 
M60A1  tank  gun  may  also  be  desirable  in  the  future  if  additional  per- 
formance specifications  are  defined  for  the  gun  or  if  conventional  tech- 
niques fail  to  achieve  the  goals  when  additional  nonlinearites  are 
included. 
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Figure  B-7.  Laminar  Vortex  Sensor  in  R 


•oooo 
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Figure  R-8.  Hull  Dynamics  in  Rate  Control  System 
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